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This  is  the  final  report  of  a three  year  project  titled,  "A  Systems  Anal  vs  if 
of  Water  Quality  Survey  Design." 

In  this  project  a study  was  made  of  water  quality  surveys  conducted  bv  the 
United  States  Army  Environmental  Hygiene  Agency  (AF.HA) . Mainly  data  and  reports 
from  studies  of  Army  Ammunition  Plants  (AAF)  were  used. 

The  focus  of  this  project  was  the  development  of  computer  aided  procedures 
which  would  assure  efficient  use  of  manpower  and  equipment  and  .ssure  that  the 
measurements  taken  give  a reasona lbe  representation  of  the  system.  l’lannin 
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survey,  conducting  the  survey  and  reporting  on  the  survey  were  included  in  the 
study. 

The  site  modeling  program  models  the  manufacturing  processes  which  contri- 
bute pollutants  to  the  system,  models  the  sewer  system,  and  models  the  treat- 
ment system  including  acid  or  caustic  neutralization,  settling  ponds,  and 
domestic  treatment.  The  inputs  to  the  model  are  the  production  levels  of  the 
manufacturing  processes  and  the  outputs  are  the  predicted  pollutant  measurement 
values  at  each  possible  measure  point  in  the  system. 

The  resource  matching  program  accepts  data  defining  proposed  measurements 
and  matches  these  against  the  available  time,  manpower,  and  equipment.  The 
output  lists  the  pollutant  to  be  measured  at  each  measure  point,  the  total 
commitment  of  time  for  each  analyst  and  for  each  piece  of  equipment.  Note  is 
made  of  any  overcommitment  of  manpower  or  equipment. 

The  model  refinement  or  updating  program  accepts  measurements  taken  durini 
a preliminary  survey  or  during  a regular  survey  and  computes  suggested  new 
parameters  for  the  process  models. 

The  indicator  model  program  evaluates  the  performance  of  sanitary  treat- 
ment facilities. 

The  program  uses  design  data,  data  from  the  operating  log  and/or  data 
generated  during  the  survey  and  computes  key  operational  characteristics.  Com 
paring  these  with  desirable  values  as  cited  in  design  books  and  manuals  will 
give  the  survey  planner  'nsight  into  the  operation  of  the  system  and  suggest 
the  need  for  more  sui  /cy  measurements  or  the  need  for  changes  in  operation. 

A system  was  developed  for  automatic  instrumentation  of  pH,  conductivity, 
and  other  parameters  which  use  strip  charp  recordings.  Interface  hardware  was 
selected  and  purchased  and  interface  software  was  developed  for  direct  connec- 
tion to  a digital  computer. 

A data  handling  syst’em  was  developed  for  use  during  and  after  the  survev. 
A PDP8-OS/8  and  peripheral  equipment  was  purchased.  Software  was  developed  to 
perform  data  handling  functions  and  to  direct  the  user  in  application  of  the 
software.  The  program  accepts  raw  data  from  the  analytical  chemist  and  per- 
forms data  conversions,  transcriptions,  and  data  logging  functions.  Output 
is  available  in  several  forms  as  may  be  needed  for  various  reports  during  and 
at  the  end  of  the  survey. 

Recommendations  are:  the  survey  planner  should  obtain  sufficient  data  in 

a preliminary  survey  to  model  and  analyze  the  site;  measurements  should  be  aut 
mated  to  the  maximum  extent  possible;  data  handling  should  be  delegated  to  t h< 
computer  when  the  operations  are  well  defined  and  repetitive.  The  programs, 
software  and  hardware  included  here  will  assist  the  survey  planner  in  followin 
these  recommendations  and  design  a more  effective  survey. 
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ABSTRACT 


The  U.S.  Army  Environmental  Hygiene  Agency's  Sanitary  Engineer  is 
tasked  with  making  decisions  in  a wide  variety  of  water  and  wastewater 
treatment  areas.  To  assist  him  in  two  of  these  areas,  this  project  was 
undertaken  with  an  objective  of  developing  some  computer  tools  with  the 
engineer  in  mind. 

Estimation  of  a wastewater  treatment  plant's  capability  to  cope  with 
an  increased  load  due  to  a post's  expansion  is  the  first  area.  Given  the 
approximate  character  of  the  wastewater  resulting  from  a post's  expansion, 
a steady-state  model  of  the  treatment  process  is  used  to  predict  the 
character  of  the  treatment  plant's  effluent.  Some  guidance  is  thus  provided  to 
to  the  engineer  in  determining  the  present  treatment  scheme's  ability  to 
provide  adequate  treatment  after  expansion. 

The  other  area  in  which  aid  is  provided  is  in  the  direction  and  analysis 
of  the  water  quality  engineering  special  study  or  ^survey*.  Utilizing 
design  and  operational  data,  a computer  model  evaluates  the  operational 
characteristics  of  the  system.  Evaluation  may  be  accomplished  with 
operational  data  available  before  the  survey  and  with  the  data  generated 
by  the  survey.  Operational  characteristics  calculated  are  then  used  by 
the  engineer  to  evaluate  the  system’s  operational  condition  and  to  then 
plan  and  conduct  the  survey. 
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wastewater  plant  serving  the  installation  will  experience  fluctuations  in 
loading.  In  addition  to  these  changes  in  personnel  strength,  many  other 
adjustments  to  the  operation  of  the  post  can  have  dramatic  effects  on  the 
loading  of  the  wastewater  treatment  plant. 

Projecting  the  effect  of  such  new  loadings  on  a wastewater  treatment 
facility,  the  sanitary  engineer  must  make  decisions  as  to  what  physical  or 
operational  changes  may  be  required  to  ensure  continued  satisfactory  plant 
performance.  Estimates  of  the  new  loading  produced  by  the  population  or 
other  changes  can  be  made  by  the  engineer  using  one  of  the  several  methods 
available  in  the  sanitary  engineering  literature.  An  estimate  of  the  future 
loadings  in  itself,  nonetheless,  does  not  establish  the  plant's  ability  to 
perform  satisfactorily  under  the  new  loading.  The  new  loading  estimate, 
however,  can  be  used  as  an  input  to  a computer  model  of  the  treatment  plant 
to  predict  the  plant's  future  performance.  Such  predictions  can  be  employed 
by  the  sanitary  engineer  in  formulating  a plan  of  action  to  meet  future  re- 
quirements. The  model  will  not  give  exact  predictions  of  the  plant  effluent. 
If  the  engineer  is  cognizant  of  this  fact,  however,  the  model  will  serve  as 
a valuable  aid  in  the  decision  making  process. 

The  TSM  simulates  two  basic  types  of  treatment  systems:  those  employ- 

ing either  (1)  a trickling  filter,  or  (2)  an  activated  sludge  system.  These 
two  types  were  chosen  because  of  their  prevalent  use  at  U.S.  Army  installa- 
tions . 


Indicator  Model 


Analysis  of  the  operation  of  wastewater  treatment  plants  is  a necessity 
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when  increased  loadings  are  expected.  Of  equal  importance  is  the  surveil- 
lance and  evaluation  of  a plant's  routine  operations.  Subsequent  recom- 
mendations regarding  process  modifications  ensure  continued  satisfactory 
performance.  Consequently,  a detailed  examination  of  the  system's  opera- 
tional condition  is  periodically  required.  This  is  accomplished  by  the 
sanitary  engineer  through  the  water  quality  engineering  special  study  or 
"survey" . 

In  preparation  for  a survey,  the  sanitary  engineer  or  "survey  officer" 
must  decide  on  the  protocol  for  the  survey,  that  is,  what  water  quality 
parameters  are  of  most  interest  and  where  in  the  system  do  they  need  to  be 
determined.  This  is,  of  course,  predicated  on  what  information  is  needed 
to  describe  the  operational  state  of  the  system.  Models  have  been  developed 
by  others  for  this  purpose.  In  particular,  there  are  the  Process,  Topology, 
and  Resource  Matching  Models  developed  by  the  Clemson  University  Department 
of  Electrical  and  Computer  Engineering.  One  feature  of  the  Topology  model 
is  the  provision  allowing  the  survey  officer  to  flag  for  measurement  during 
the  survey,  any  particular  water  quality  parameter  at  any  point  in  the  sys- 
tem. Flaging  in  this  manner  usually  is  done  when  the  survey  officer  "feels" 
a particular  parameter  may  have  a significant  effect  upon  the  overall  opera- 
tional state  of  the  system. 

A preliminary  evaluation  of  the  process  variables  of  a treatment  plant 
will  provide  an  initial  estimate  of  its  operational  state.  These  process 
variables,  such  as  weir  overflow  rates  and  process  loading  intensity,  are 
functions  of  the  design  and  operation  of  the  treatment  plant  and  the  char- 
acter of  the  wastewater.  Thus,  a preliminary  evaluation  of  the  process 
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variables  will  indicate  possible  areas  where  stress  on  the  plant  exists. 

The  survey  officer  then  determines  whether  certain  key  parameters  should  be 
flaged  and  whether  measurements  within  the  wastewater  treatment  plant  should 
be  made  during  the  survey. 

After  completion  of  the  survey,  the  survey  officer  must  analyze  the 
data  collected.  This  analysis  results  in  the  final  determination  of  the 
treatment  plant's  operational  condition  and  forms  a basis  for  the  survey 
officer's  recommendations  to  improve  treatment  plant  performance.  Again, 
evaluation  of  the  process  variables,  this  time  using  the  data  generated  by 
the  survey,  will  assist  the  survey  officer  in  determining  the  plant's  opera- 
tional condition.  These  variables  will  indicate  areas  where  design  or  opera- 
tional inadequacies  exist  which  may  be  the  cause  of  unsatisfactory  plant 
performance.  Trouble  spots  in  the  system  can,  therefore,  be  pinpointed. 

Also,  the  survey  officer  can  determine,  through  analysis  of  the  process 
variables,  what  remedial  action  will  likely  improve  plant  performance. 

The  Indicator  Model  was  developed  to  evaluate  these  process  variables 
for  two  types  of  wastewater  treatment  systems:  trickling  filter  and  activated 
sludge  type  of  plants.  The  model  can  be  used  to  compute  the  process  variables 
utilizing  preliminary  operational  data  that  is  available  prior  to  the  incep- 
tion of  the  survey  or  the  operations  data  generated  during  the  survey. 
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II  TREATMENT  SYSTEM  MODEL 
Introduction 

A wide  variety  of  computer  models  have  been  developed  to  simulate  dif- 
ferent wastewater  treatment  systems.  Both  steady  state  and  dynamic  models 
have  been  proposed.  While  the  steady  state  model  does  not  represent  the 
transient  operation  of  a process  and,  as  such,  is  generally  more  simplistic 
than  the  dynamic  model,  it  does  give  an  estimate  of  a process'  long-term 
performance.  Since  this  long-term  estimate  is  what  the  model  is  intended  to 
provide,  the  steady  state  approach  has  been  used  in  this  development.  To 
quote  Andrews  [12]: 

...the  simplest  possible  expression  should  be  used.  The  pur- 
poses for  which  the  model  is  to  be  used  must  be  defined  so  that 
one  can  be  developed  which  is  adequate  for  the  intended  use. 

Mathematical  elegance  should  be  secondary  and  it  should  be  re- 
membered that  a model  which  is  too  complex  may  be  subject  to  either 
misuse  or  disuse. 

Model  Formulation 

Actually,  the  overall  model  which  has  been  developed  is  a combination  of 
several  steady  state  process  models.  All  the  unit  processes  that  are  includ- 
ed in  the  system  have  been  modeled  and  linked  one  to  another.  These  include: 
primary  clarification,  secondary  clarification,  the  trickling  filter,  and  the 
activated  sludge  process.  As  mentioned  before,  these  unit  processes  are 
grouped  into  two  basic  treatment  trains,  one  incorporating  the  trickling 
filter  while  the  other  includes  the  activated  sludge  process  (see  Figure  1). 
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a.  TRICKLING  FILTER  TRAIN 


b.  ACTIVATED  SLUDGE  TRAIN 
Figure  1:  BASIC  TREATMENT  TRAINS 


Complex  interactions  between  unit  processes  are  not  explicitly  considered  in 
the  formulation  of  the  model.  Again,  this  is  not  necessary  for  the  intended 
use  of  this  treatment  system  model.  A more  detailed  explanation  of  these  com- 
plex interactions  may  be  found  in  the  sections  dealing  with  the  two  treatment 
trains. 


Water  Quality  Evaluation 


Trickling  filter  and  activated  sludge  systems  are  both  biological  pro- 
cesses, i.e.  , they  employ  microorganisms  to  assimilate  and  metabolize  the  pol- 
lutants in  the  wastewater.  As  such,  these  two  systems  can  only  be  used  to  treat 
biodegradable  wastes  su~-i  as  "domestic"  wastes  and  certain  types  of  industrial 
wastes.  To  establish  the  quality  of  these  wastes,  a number  of  water  quality 
parameters  have  been  used  in  the  past. 

Four  such  parameters  are  most  commonly  used  by  sanitary  engineers,  not 
only  to  indicate  the  wastewater  quality,  but  also  to  determine  the  amount  of 
treatment  accorded  the  wastewater  by  the  treatment  system.  These  parameters 
are  as  follows:  1)  five-day  biochemical  oxygen  demand  (B0D5);  2)  suspended 

solids  (SS);  3)  ammonia;  and  4)  nitrates.  Given  the  influent  levels  of  these 

four  water  quality  parameters  and  the  influent  flow  rate,  the  model  predicts 
their  respective  concentrations  in  the  effluent  of  a treatment  plant. 

Performance  of  each  unit  process  which  is  incorporated  in  the  treatment 
train  is  described  by  a mathematical  expression.  The  expression  is  used  to 
calculate  the  values  of  the  water  quality  parameters  in  the  effluent  of  the 
process.  These  values  then  become  the  input  to  the  next  unit  process  in  the 
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treatment  train.  The  output  conditions  calculated  for  the  last  unit  process 


in  the  treatment  train  become  the  plant's  effluent  conditions. 

Effects  of  Process  Dynamics 

Since  the  TSM  is  a steady  state  model,  the  influent  and  predicted  effluent 
concentrations  must,  of  necessity,  be  averages.  If  the  influent  flow  rate  or 
concentrations  are  known  to  fluctuate  drastically,  the  steady  state  model  pre- 
dictions are  less  accurate.  The  engineer  must  consider  this  fact  when  apply- 
ing the  model  and  when  interpreting  the  results  of  the  model. 

Trickling  Filter  System 

Was,.:water  treatment  systems  utilizing  the  trickling  filter  can  take  on 
many  different  configurations  depending  on  whether  the  filters  are  staged  and 
on  the  recycle  pattern.  The  model  presented  here  allows  for  three  distinct 
variations.  All  three  are  for  single-stage  filters,  i.e.  , filters  in  a parallel 
but  not  in  series.  The  first  is  for  no  recycle;  the  second  for  recycle  from 
the  filter  effluent;  and  the  third  for  recycle  from  the  secondary  clarifier 
effluent  (see  Figures  2 a,b,c). 

The  type  of  system  modeled  is  controlled  by  specification  of  a recircula- 
tion data  input  to  the  model.  If  both  of  the  recycle  terms  are  initialized 
at  zero,  then  a no  recycle  system  is  specified.  (Figure  2a).  If  the  recycle 
to  the  filter  is  not  zero,  but  recycle  from  the  filter  effluent  is,  then  a 
system  with  recycle  from  the  secondary  clarifier  will  be  modeled  (Figure  2c). 
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. Trickling  Filter:  No  Recycle 


. Trickling  Filter:  Recycle  from  Trickling  Filter  Effluent 


c.  Trickling  Filter:  Recycle  from  Secondary  Clarifier  Effluent 


d.  Activated  Sludge  System 


Figure  2:  Schematic  of  Treatment  Systems 


Finally,  if  the  recycle  flow  rate  from  the  filter  is  equal  to  the  recycle 
flow  rate  to  the  filter  and  is  not  equal  to  zero,  the  model  simulates  a system 
where  recycle  is  from  the  filter  effluent  (Figure  2b). 

Other  Considerations 

Factors  which  may  effect  the  operation  of  the  system  but  are  not  incor- 
porated in  the  formulation  of  the  model  include  the  effect  of  the  trickling 
filter  media  type  on  the  BOD  reduction;  certain  complex  interactions  between 
the  unit  processes;  and  the  additional  hydraulic  and  solids  loading  imposed 
on  the  primary  clarifier  by  wasting  secondary  sludge  to  the  primary  clarifier 
influent. 

One  example  of  complex  interactions  between  unit  processes  is  the  effect 
of  the  primary  clarifier  on  the  performance  of  the  trickling  filter.  The  model 
predicts  the  case  in  which  the  primary  clarifier  is  overloaded  thereby  per- 
mitting large  suspended  solids  carryover.  Conversely,  the  model  does  not  ac- 
count for  the  situation  if  the  high  solids  concentration  in  the  trickling 
filter  influent  were  to  clog  the  filter  or  otherwise  degrade  its  performance. 
While  the  exclusion  of  these  and  other  considerations  may  result  in  minor 
errors,  such  errors  are  usually  negligible.  With  an  understanding  of  the 
structure  of  the  model,  as  just  presented,  and  a knowledge  of  the  mathematical 
relationships  incorporated  in  the  model,  the  engineer  can  utilize  the  present 
trickling  filter  system  model  to  great  advantage  or  modify  it  to  more  suitably 
fit  his  needs. 
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Mathematical  Relationships 


A listing  of  the  expressions  that  were  used  in  forming  the  overall  model 
of  the  trickling  filter  system  follows.  These  relationships  are  used  to  eval- 
uate the  operation  of  the  unit  processes. 


Primary  Clarifier 


RF  = 1.0  - 0.82EXP(-Q/2780A) 
1.0-  0. 00205EXP( -Q/2780A) 


RF  = Removal  factor 
Q = Plant  flow  (MGD) 

A = Total  surface  area  of  clarifier  (acres) 


SS  = SS.  * RF 
out  in 


B0Dout  = (BODin)(0.3)  +(B0Din)(RF)(0.7) 


No  effect  on  ammonia  or  nitrates 


Trickling  Filter 


BOD  reduction  L J 

E = 1.0 

1.0+0.0085/  WU 


E = Decimal  BOD  removal  efficiency 
W = BOD  loading  (lbs/day)  to  filter 
(not  including  recycle) 

V = Volume  of  filter  (Acre  - Ft) 

F = Recirculation  factor  1 .0+R  where  R = recycle 

( 1 . 0+0 . 1 R)  2 influent 


i 

! 


! 


A temperature  correction  has  also  been 
incorporated  [3],  [4] 

T-20 

Et  = E (1 .038) 


E.  = Filter  efficiency  at  temperature,  T 
El  = Filter  efficiency  as  above 
T = Temperature  of  wastewater,  deg  C 


E°°out  = B001n*  (’.0  - Et) 
Suspended  Solids 


Lbs  solids  produced  = 0.2(lbs  BOD.  -lbs  BOD  .) 

in  out7 

SS  . = SS.  + lbs  solids  produced 
8 . 34 ( f 1 ow+recyc  1 e ) 

Nitrification^ 

E = EXP(  -KD/Qn) 

E = Decimal  reduction  in  NH^-N  concentration 
K = Reaction  constant  related  to  the  specific 
surface  of  the  media 
D = Depth  of  filter  (Ft) 

Q = Hydraulic  loading  rate  (MGAD) 
n = Constant  related  to  the  specific  surface  and 
configuration  of  packing 

K?f=K  at  20°  C =0.23  1 " rock  media 
=0.13  2i"  rock  media 

K=K20(1 ,07)T"20  where  T is  temp.  (°C) 

n = 2.36  1 " rock  media 

3.80  2i"rock  media 

[NH3-N]out  = [NH3-N]jn*E 


[N°3-»]out  ' [N03-N]jn  * [NH3-N],n  - [NH3-N]out 
Secondary  Clarifier^  ^ 

The  secondary  clarifier  for  the  trickling  filter 
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system  is  assumed  to  operate  similar  to  the 
primary  and  the  same  equation  is  utilized. 

The  exception  being  that  only  suspended  solids 
removal  is  considered. 

Data  Input 

Following  is  a listing  of  the  data  input  necessary  to  implement  the 
simulation  of  a trickling  filter  plant. 

Card  1 

Column  1:  Type  of  plant;  1 for  trickling  filter  plant 

Columns  2-9:  Flow  rate  to  plant  (MGD) 

Columns  10-17:  BOD  concentration  in  wastewater  (mg/1) 

Columns  18-25:  Suspended  solids  cone,  in  water  (mg/1) 

Columns  26-33:  Ammonia  cone,  in  water  (mg/1  NH^  as  N) 

Columns  34-41:  Nitrate  cone,  in  water  (mg/1  N03  as  N) 

Columns  42-49:  Temperature  of  water  (°C) 

Card  2 

Columns  1-8:  Total  surface  area  of  primary  clarifier  (acres) 

(2  clarifiers  @ 2 acres  ea.  = 4 acres) 

Card  3 

Columns  1-8:  Total  filter  volume  (acre-ft) 

(2  filters  0 3 acre-ft  ea.  = 6 acre-ft) 

Columns  9-16:  Total  filter  surface  area  (acres) 

(3  filters  @ 4 acres  ea.  = 12  acres) 

Columns  17-24:  Total  recycle  to  filters  (MGD) 

Columns  25-32:  Total  recycle  from  filter  effluent  (MGD) 

Media  coefficient  (K)  0.23  1"  rock  media 

0.13  2i"  rock  media 
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Columns  33-40 


Columns  41-48:  Media  coefficient  (n)  2.36  1"  rock  media 

3.80  2i"  rock  media 


Card  4 

Columns  1-8:  Total  surface  area  of  secondary  clarifier  (acres) 

(2  clarifiers  @ 2 acres  = 4.0) 


As  with  the  trickling  filter  type  of  wastewater  treatment  plant,  many  dif- 
ferent modifications  of  the  activated  sludge  process  exist.  The  mode  of  op- 
eration of  the  treatment  plant  determines  which  process  modification  it  employs. 
Of  first  concern  is  the  fluid  regime,  i.e.  , whether  the  flow  is  complete-mix 
or  plug-flow.  The  model  presented  here  is  based  on  complete-mix  conditions. 

It  can  also  be  applied  to  a plug-flow  operation  though,  since  the  complete-mix 
formulation  generally  underestimates  the  treatment  level  attainable  with  a 
plug-flow  system.  Because  this  underestimation  is  not  generally  severe,  the 
model  can  be  considered  adequate  for  application  to  plug-flow  systems.  The 
model,  therefore,  can  simulate  the  complete-mix  process  including  the  extended 
aeration  modification  of  that  process,  and  is  acceptable  for  the  plug-flow 
modification  (See  Figure  2d). 

Specification  by  the  engineer  of  the  type  of  activated  sludge  system  to  be 
simulated  by  the  model  is  not  necessary.  Input  of  the  correct  aeration  basin 
operational  characteristics  will  result  in  a simulation  of  the  appropriate 
modi f ication . 


Aeration  Basin  Operational  Characteristics 
Accuracy  of  the  model's  predictions  is  also  highly  dependent  upon  the 
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estimate  of  the  aeration  basin's  operational  characteristics  made  by  the 
survey  engineer.  To  obtain  reliable  results,  values  for  the  mixed  liquor 


suspended  solids,  recycle  rate,  and  net  yield  terms  must  closely  approximate 
those  that  exist  during  actual  operation. 

While  accurate  estimates  of  these  characteristics  are  mandatory,  it  is 
often  quite  difficult  to  obtain  them.  The  aeration  basin's  interdependence 
on  the  operation  of  the  other  unit  processes  which  surround  it,  especially  the 
secondary  clarifier,  is  the  cause  of  this  difficulty.  In  the  situation  where 
the  secondary  clarifier  is  underloaded  (that  is,  the  sludge  blanket  is  at  the 
bottom),  the  MLSS  is  determined  solely  by  the  amount  of  sludge  wasted  and  is 
easily  controlled  at  the  desired  level.  However,  future  conditions  may  result 
in  a critically  loaded  or  overloaded  clarifier.  This  situation  is  a function 
of  the  secondary  clarifier  solids  loading  and  the  settleabil ity  of  the  sludge. 

Here,  the  interdependence  between  aeration  basin  and  secondary  clarifier  exists 
and  governs  what  the  aeration  basin's  operational  characteristics  will  be. 

Should  overloading  occur,  plant  failure  will  result.  An  overload  can  not 
be  predicted  by  the  model  but  must  be  accounted  for  by  the  engineer  when  he 
estimates  the  aeration  basin's  operational  characteristics. 

One  of  the  most  useful  tools  for  estimating  these  parameters  with  ac- 
curacy is  experience.  Comparison  of  the  situation  to  those  currently  in 
existence  or  in  the  literature  that  are  similar  in  nature  can  also  be  of  con- 
siderable aid. 

With  an  understanding  of  the  structure  of  the  activated  sludge  model 
as  just  presented,  and  a knowledge  of  the  mathematical  relationships  that 

describe  the  operation  of  the  different  unit  processes  in  the  system,  the  j 
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engineer  should  be  able  to  use  the  model  to  advantage  or  modify  it  to  more 
suitably  fit  his  needs. 


Mathematical  Relationships 


Following  is  a listing  of  the  expressions  that  were  used  in  forming 


the  overall  model  of  the  activated  sludge  system. 


Primary  Clarifier 

The  relationship  is  the  same  as  that  used  for  the  trickling 
filter  system. 

Aeration  Basin 

BOD  reduction^ 


BOD 


out 


BODin 


(FlCWR) 


MLSS  = Mixed  liquor  suspended  solids  (mg/1) 
Vol  = Aeration  basin  volume  (MG) 

Yn  = Net  yield  (lb  SS  produced/lb  BOD  used) 

FLOWR  = Plant  flow  + recycle  flow  (MGD) 

Theta  = Sludge  age  or  sludge  retention  time 
= 1 .0  where  PL  I = BOD.  (Flow) 

(MLSS) ( Vol ) 

Nitrification^ 


A double  "halving"  or  Balzano  type  iteration  is  performed 
on  the  effluent  concentrations  of  ammonia  and  the  nitrify- 
ing bacteria,  Nitosomonas,  until  the  governing  system 
equations  are  satisfied: 


Si  = 


KrSoXoDt 


Xi 


K S X D 
r o o 

K$  + S 


t 

o 
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where: 


Si  = Amnonia  in 

S = Ammonia  out 
o 

Kr  = Synthesis  rate  constant 
Y = Yield  coefficient 
K$  = Saturation  coefficient 
Xi  = Nitrosomonas  in 
XQ  = Nitrosomonas  out 

Detention  time  (Vol/Flow) 

A detailed  explanation  and  justification  of  assumptions  may  be 
found  in  8. 

roi 

Secondary  Clarifier1-  J 

SS  = 4.5  + 8.6SSF 
out 

SSF  = SSin  (FR)  3.785  = Solids  surface  feed  (kg/m2  - hr) 
24A 

where:  FR  = Plant  flow  + recycle  (MGD)  „ 

A = Surface  area  of  clarifier  (nr) 


Data  Input 


Following  is  a listing  of  the  data  input  necessary  to  implement  the 
simulation  of  an  activated  sludge  plant. 

Card  1 

Column  1:  Type  of  plant:  2 for  activated  sludge  plant 

Columns  2-9:  Predicted  flow  to  plant  (MGD) 

Columns  10-17:  Predicted  BOD  of  influent  (mg/1) 

Columns  18-25:  Predicted  suspended  solids  of  influent  (mg/1) 
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Columns  26-33:  Predicted  ammonia  in  influent  (mg/1  NH^  as  N) 
Columns  34-41:  Predicted  nitrate  in  influent  (mg/1  NO^  as  N) 
Columns  42-49:  Predicted  temperature  of  influent  (deg.  C) 

Card  2 

Columns  1-8:  Total  surface  area  of  primary  clarifiers  (acres) 

(2  clarifiers  @ 2 acres  ea.  = 4 acres) 

Card  3 

Columns  1-8:  Total  volume  of  aeration  tanks  (MG) 

2 tanks  @ 3 MG  ea.  = 6 MG 

Columns  9-16:  Estimated  operating  concentration  of 

MLSS  (mg/1) 

Columns  17-24:  Estimated  operating  recycle  rate  (MGD) 

Columns  25-32:  Net  yield  (lb  SS/lb  BOD) 

If  not  readily  available,  calculate: 

Y = (1b  SS  wasted/day)  + (effluent  SS)  (Flow) 
tB0D,„  - ““out1  F,0“ 

Card  4 

Columns  1-8:  Total  surface  area  of  secondary  clarifiers  (acres) 

(2  clarifiers  @ 2 acres  ea.  = 4 acres) 
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Ill  INDICATOR  MODEL 


Introduction 


The  indicator  model  is  used  to  calculate  certain  parameters,  gen- 
erally termed  process  variable  or  operating  characteristics,  that  indicate 
the  operational  state  of  the  system.  These  calculated  values  can  then  be 
compared  with  the  desirable  values  as  cited  in  design  books  or  manuals 
[7,  10,  11].  A comparison  will  give  the  engineer  insight  into  the  oper- 
ation of  the  system  such  that  problem  areas  and  possible  causes  of  these 
problems  are  more  readily  apparent.  It  can  also  indicate  whether 
inadequate  plant  performance  is  due  to  design  deficiencies  or  poor 
operational  practices. 

An  example  should  help  to  illustrate  the  use  of  the  indicator  model 
more  clearly.  If  a plant's  effluent  suspended  solids  concentration  is  too 
high,  the  engineer  might  first  suspect  the  secondary  clarifier.  Solids 
removal  efficiency  calculated  by  the  model  will  indicate  whether  or  not 
the  secondary  clarifier  is  performing  poorly.  If  the  secondary  clarifier 
does  appear  to  be  performing  poorly,  operating  characteristics  such  as 
surface  settling  rate  and  solids  loading  rate,  calculated  by  the  model, 
can  be  compared  to  those  suggested  in  the  literature.  If  the  clarifier 
is  failing  due  to  a hydraulic  overload,  comparison  of  the  surface  settling 
rate  to  common  values  should  indicate  this.  If  the  operational  character- 
istics of  the  secondary  clarifier  compare  favorably  with  those  in  the 
literature,  but  it  still  appears  to  be  failing,  an  investigation  of  the 
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operational  characteristics  of  the  other  unit  processes  might  be  helpful. 
Should  the  process  loading  intensity  of  the  aeration  basin  not  compare 
favorably  with  cited  values,  a poorly  settling  sludge  might  well  be  the 


cause  of  the  secondary  clarifier's  failure. 

I 

Operational  Characteristics 

Like  the  treatment  system  model,  the  indicator  model  is  applicable 
to  trickling  filter  and  activated  sludge  type  treatment  plants.  Opera- 
tional characteristics  calculated  by  the  model  are  as  follows: 

Primary  Clarifier 
Detention  time  (hrs) 

Surface  settling  rate  (GPD/ft^) 

Weir  overflow  rate  (GPD/ft) 

Solids  loading  rate  (lb  SS/ft^-day) 

Solids  removal  efficiency  {%) 

Trickling  Filter 
Hydraulic  loading  (MGAD) 

Organic  loading  (lb  BOD/day/acre-ft) 

Recirculation  ratio 
BOD  removal  efficiency  (%) 

Aeration  Basin 
Detention  time  (hrs) 

Volumetric  loading  (lb  B0D/day/1000  ft'*) 

Process  loading  intensity  (lb  BOD/day/lb  MLSS) 
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Sludge  age  or  retention  time  (days) 

Diffused  aeration 

Volumetric  air  supply  (CFM/1000  ft3) 
Organic  air  supply  (ft3/lb  BOD) 
Mechanical  aeration 

Volumetric  power  supply  (HP/1000  ft3) 
Organic  oxygen  supply  (lb  oxygen/lb  BOD) 
BOD  removal  efficiency  (%) 

Secondary  Clarifier 
Detention  time  (hrs) 

Surface  settling  rate  (GPD/ft3) 

Weir  overflow  rate  (GPD/ft) 

Solids  loading  rate  (lb  SS/day/ft3) 

Solids  removal  efficiency  (%) 

Overall  Plant 

BOD  removal  efficiency  (%) 

Solids  removal  efficiency  (%) 


Mathematical  Relationships 


Each  of  the  operational  characteristics  calculated  by  the  model  is 
represented  by  a mathematical  expression.  Many  of  these  expressions  are 
evaluated  at  both  average  and  peak  flow  conditions.  This  is  done  by 
using  either  the  average  or  peak  flow  term  in  the  expression.  However, 
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this  may  not  always  be  completely  accurate,  as  the  BOD  suspended  solids 
concentrations  also  vary  diurnally.  Peak  values  for  the  operational 
characteristics,  therefore,  may  be  in  error  because  only  the  peak  flows 
and  not  the  peak  concentrations  were  used  in  the  calculations.  Such 
error  may  be  overcome,  if  desired,  by  substituting  the  peak  values  of 
these  parameters  (BOD,  flow,  etc.)  for  the  average  values  as  input  data 
and  running  the  model  another  time.  Values  of  the  operational  character- 
istics calculated  on  this  second  run,  which  are  printed  out  as  average, 
will  then  represent  the  peak  values.  Following  is  a listing  of  the 
relationships : 

Primary  Clarifier 

Detention  time:  DT  = V 

1000  F/24 

V = clarifier  volume  (1000  gal) 

F = plant  flow  (MGD) 

Surface  settling  rate: 

SSR  = F * 106 

B 

F = plant  flow  (MGD)  2 

A = clarifier  surface  area  (ft  ; 

Weir  overflow  rate: 

WOR  = F * 106 

E 

F = plant  flow  (MGD) 

L = weir  length  (ft) 

Solids  loading: 

SL  = (SS)  (F)  (8.34) 

SS  = influent  suspended  solids  (mg/1) 

F = plant  flow  (MGD)  2 

A = clarifier  surface  area  (ft  ) 
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Solids  removal  efficiency: 


1 


SRE  = SSI  - SSO 

— 35T — 


* 100. 


SSI  = plant  Influent  suspended  solids  (mg/1) 
SSO  = primary  effluent  suspended  solids  (mg/1) 


Trickling  Filter 
Hydraulic  loading: 

HL  = F + R 

~T~ 

F = plant  flow  (MGD) 

R = recycle  flow  to  filter  (MGD) 
A = filter  surface  area  (acres) 

Organic  loading: 

OL  = (F)  (BOD)  (8.34) 

V 

F = plant  flow  (MGD) 

BOD  = primary  effluent  BOD  (mg/1) 

V = filter  volume  (acre-ft) 

Recirculation  ratio: 

RR  = R 

T 

F = plant  flow  (MGD) 

R = recycle  flow  to  filter  (MGD) 

BOD  removal  efficiency: 

E = B0D1  - BODq  * 100> 

TOITj 

B0D,=  primary  effluent  BOD  (mg/1) 
B0Dq=  plant  effluent  BOD  (mg/1) 

Aeration  Basin 
Detention  time: 

DT  = 24V 

(F+ft)TO 
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A 


Aeration  Basin  (Continued) 


V = aeration  basin  volume  (1000) 

F = plant  flow  (MGD) 

R = recycle  flow  (MGD) 

Volumetric  loading: 

VL  = (BOD)  (F)  (8.34)  (7.481) 

V 

BOD  = primary  effluent  BOD  (mg/1) 

Process  loading  intensity: 

pu-«ifr 

BOD  = primary  effluent  BOD  (mg/1) 

F = plant  flow  (MGD) 

MLSS=  mixed  liquor  suspended  solids  (mg/1) 

V = aeration  basin  volume  (1000  gal) 

Sludge  age: 


(MLSSH 

0.001V) (8. 34) 

(F  }(RSS)(8.34 

) + ( ) (Ss) ( &. 34) 

MLSS=  mixed  liquor  suspended  solids  (mg/1) 

V = aeration  basin  volume  (1000  gal) 

Fw  = waste  sludge  flow  rate  (MGD) 

RSS  = return  sludge  suspended  solids  (mg/1) 

F = plant  flow  (MGD) 

SS  = plant  effluent  suspended  solids  (mg/1) 
A basic  assumption  made  in  the  sludge  age  formulation  is  that 
sludge  wasting  is  from  the  sludge  recycle  line.  If  the  wasting 
is  from  the  aeration  basin,  the  actual  waste  sludge  flow  rate 
must  be  adjusted  before  input  to  the  model. 

F,  • (F.ymss) 

F.  = input  waste  sludge  flow  (GPD) 

F,  = actual  waste  sludge  flow  (GPD) 

d 

MLSS=  mixed  liquor  suspended  solids  (mg/1) 
RSS  = return  sludge  suspended  solids  (mg/1) 
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Volumetric  air  supply: 


f 


VAS 


! 

1 


AV 

V 


Organic  air  supply: 


OAS 


AV 

BOD 

F 


( AV) ( 1000) (7.481) 

immy) 

3 

volume  of  air  delivered  (1000  ft  /day) 
aeration  basin  volume  (1000  gal) 


3 

volume  of  air  delivered  (1000  ft  /day) 
primary  effluent  BOD  (m^/1) 
plant  flow  (MGD) 


Volumetric  power  supply: 

VPS  = 7.481  HP 
V 

HP  = total  aerator  horsepower  (HP) 

V = aeration  basin  volume  (1000  gal) 

Organic  oxygen  supply: 


N = oxygen  transfer  rate  (lb  02/HP-hr) 
HP  = total  aerator  horsepower  (HP) 

HR  = hours  of  aerator  operation  (hr/day) 
BOD  = primary  effluent  BOD  (mg/1) 

F = plant  flow  (MGD) 

N = (2.5)  CS-C  (1.024)  T-20  (0.8) 

07 


CS  = D.O.  saturation  in  aeration  basin  (mg/1) 
C = D.O.  cone,  in  aeration  basin  (mg/1) 

T = plant  effluent  temperature  (°C) 

A basic  assumption  made  in  the  organic  oxygen  supply  formulatior  is 

that  the  mechanical  aerators  employed  are  rated  at  2.5  lb  02/HP-hr 

under  standard  conditions  as  defined  in  [7].  Also  assumed  is  an  alpha 
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value  for  the  wastewater  of  0.8.  If  the  engineer  has  the  actual  values 
of  these  parameters  for  the  system  being  investigated,  the  model  should 
be  changed  accordingly. 

BOD  removal  efficiency: 

E ■ BOD,  - BOD0  . 

^*0 

BOD..  = primary  effluent  BOD  (mg/1) 

B0Dq  = plant  effluent  BOD  (mg/1) 

Secondary  Clarifier 

For  the  secondary  clarifier,  the  flow  to  the  clarifier  depends  on 
whether  it  is  part  of  a trickling  filter  or  activated  sludge  type 
system.  For  the  trickling  filter  system,  the  flow  to  the  clarifier 
is  equal  to  the  plant  flow  plus  the  recycle  to  the  filter  minus  the 
recycle  from  the  filter  effluent.  For  the  activated  sludge  system, 
the  flow  to  the  clarifier  is  equal  to  the  plant  flow  plus  the  re- 
cycle sludge  flow. 

Detention  time: 

DT  = V 

1000  F/24 

V = clarifier  volume  (1000  gal) 

F = flow  to  clarifier  (MGD) 

Surface  settling  rate: 

SSR  = F * 106 

— a 

F = flow  to  clarifier  (MGD)  - 
A = clarifier  surface  area  (ft  j 


L _ 


Weir  overflow  rate: 


WOR  = F * 106 
[ 

F = flow  to  clarifier  (MGD) 

L = weir  length  (ft) 

Solids  loading  rate: 

SL  = (SS) (F) (8.34) 

A 

SS  = aeration  basin  or  trickling  filter 
effluent  suspended  solids  (mg/1) 

F = flow  to  clarifier  (MGD) 

Solids  removal  efficiency: 

E - SSI  - SSO  * ,nn 
— SST5 — 100 

SSI  = aeration  basin  or  trickling  filter 
effluent  suspended  solids  (mg/1) 

SSO  = plant  effluent  suspended  solids  (mg/1) 


Data  Input 

The  indicator  model  is  designed  to  use  data  that  is  available  prior 
to  a survey's  implementation  and  also  the  data  generated  by  the  survey. 
Design  data  needed  should,  of  course,  always  be  available.  Operational 
data  required  will  be  available  from  the  daily  and  monthly  operation  logs 
as  required  by  Technical  Manual  5-665,  and  then  from  the  data  generated 
by  the  survey.  (If  the  value  of  any  parameter  is  not  available,  it 
should  be  replaced  by  - 1.0  in  the  appropriate  columns  for  the  data  input. ) 
A listing  of  the  data  input  necessary  to  implement  the  indicator  model 
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follows: 


Trickling  Filter  Type  Plant 
Card  1 


Column  1:  Plant  type;  1 for  trickling  filter  type  plant. 

Columns  2-41:  Name  of  installation  e.g.  Fort  Ord,  CA 

(Should  be  centered) 


Card  2 - Primary  Clarifier 


Columns  1-8:  Total  volume  (1000  gal) 

(2  clarifiers  0 4000  gal  ea.  = 8.0) 

2 

Columns  9-16:  Total  surface  area  (Ft„) 

(2  clarifiers  @ 200  Ft^  ea.  = 400.) 

Columns  17-24:  Total  weir  length  (Ft) 

(2  clarifiers  @ 500  Ft  ea.  = 1000.) 


Columns  25-32:  Design  flow  rate  (MGD) 

Columns  33-40:  Design  detention  time  (Hrs) 

If  not  readily  available,  calculate 

DT  = Total  Vol  (1000  gal)  X 0.024 
Design  flow  (MGD) 

p 

Columns  41-48:  Design  surface  settling  rate  (GPD/Ft  ) 
If  not  readily  available,  calculate 

SSR  = Design  flow  (MGD)  X 106q 
Total  surface  area  (Ftc ) 


Columns  49-56:  Design  weir  overflow  rate  (GPD/Ft) 
If  not  readily  available,  calculate 

WOR  = Design  flow  (MGD)  X 106 
Total  weir  length  (Ft) 


Card  3 - Trickling  Filter 

Columns  1-8:  Design  flow  rate  (MGD) 

Columns  9-16:  Design  organic  loading  (Lb  BOD/Acre-Ft/Day) 
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Columns  17-24:  Total  volume  (Acre-Ft) 

(2  filters  0 3 acre-ft  ea.  = 6.) 

Columns  25-32:  Total  area  (Acres) 

(2  filters  @ .75  acres  ea.  = 1.5) 

Columns  33-40:  Total  design  recycle  per  day  (MG) 

If  not  readily  available,  calculate 

R = Design  percent  recycle  X 100.  X Design  flow  (MGD) 


Card  4 - Secondary  Clarifier 

Columns  1-8:  Total  volume  (1000  Gal) 

(2  clarifiers  @ 4000  gal  ea.  = 8.0) 

p 

Columns  9-16:  Total  surface  area  (Ft?) 

(2  clarifiers  @ 200  Ftd  ea.  = 400.) 


Columns  17-24:  Total  weir  length  (Ft) 

(2  clarifiers  0 500  Ft  ea.  = 1000.) 

Columns  25-32:  Design  flow  rate  (MGD) 

Columns  33-40:  Design  detention  time  (Hrs) 

If  not  readily  available,  calculate 


0T  = Total  Vol  (1000  gal)  X 0.024 
Design  flow  (M5D) 

2 

Columns  41-48:  Design  surface  settling  rate  (GPD/Ft  ) 
If  not  readily  available,  calculate 

SSR  = Design  flow  (MGD)  X lpjj 
Total  surface  area  (?tc) 


Columns  49-56:  Design  weir  overflow  rate  (GPD/Ft) 

If  not  readily  available,  calculate 

W0R  = Design  flow  (MGD)  X 106 
total  weir  length  (Ft) 

2 

Columns  57-64:  Design  solids  loading  (Lb  SS/Ft  -Day) 


Card  5 - General  Operations 

Columns  1-8:  Peak  flow  (influent)  (MGD) 

Columns  9-16:  Total  daily  flow  (influent)  (MG) 
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Columns  17-24:  BOD  (influent)  (mg/1) 

Columns  25-32:  Suspended  solids  (influent)  (mg/1) 

Columns  33-40:  BOD  (effluent)  (mg/1) 

Columns  41-48:  Suspended  solids  (effluent)  (mg/1) 

Card  6 - Process  Operations 

Columns  1-8:  Primary  clarifier  effluent  BOD  (mg/1) 

Columns  9-16:  Primary  effluent  suspended  solids  (mg/1) 

Columns  17-24:  Trickling  filter  effluent  BOD  (mg/1) 

Columns  25-32:  Filter  effluent  suspended  solids  (mg/1) 

Columns  33-40:  Total  recycle  to  filter  (MGD) 

Columns  41-48:  Total  recycle  from  filter  effluent  (MGD) 
Columns  49-56:  Recycle  from  secondary  clarifier  (MGD) 
Activated  Sludge  Type  Plant 
Card  1 

Column  1:  Plant  type;  2 for  activated  sludge  type  plant. 

Columns  2-41:  Same  as  for  trickling  filter  type  plant. 

Card  2 - Primary  Clarifier 

Same  as  for  trickling  filter  type  plant. 

Card  3 - Aeration  Tanks 

Columns  1-8:  Total  tank  volume  (1000  Gal) 

(2  tanks  @ 50,000  gal  ea.  = 100.) 

Columns  9-16:  Design  flow  rate  (MGD) 

Columns  17-24:  Design  detention  time  (Hrs) 

If  not  readily  available,  calculate 

DT  = Total  Vol  (1000  Gal)  X 0.024 
besign  flow  (MgD) 
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Columns  25-32:  Design  organic  loading  (Lb  B0D/day/1000  Ft^) 
Columns  33-40:  Design  mixed  liguor  suspended  solids  (mg/1) 


Columns  41-48:  Number  of  aerators  or  compressors 

Columns  49-56:  Total  aerator  or  compressor  horsepower 
(2  aerators  @ 50  HP  ea.  = 100.) 

Card  4 - Secondary  Clarifier 

Same  as  for  trickling  filter  type  plant. 

Card  5 - General  Operations 

Same  as  for  trickling  filter  type  plant  with  the  exception  of: 
Columns  49-56:  Average  temperature  of  plant  effluent  ( C) 
Columns  57-64:  D.O.  saturation  in  aeration  basin  (mg/1) 

Card  6 - Process  Operations 
Columns  1-8: 

Columns  9-16: 

Columns  17-24: 

Columns  25-32: 

Columns  33-40: 


Columns  41-48: 
Columns  49-56: 
Columns  57-64: 
Columns  65-72: 
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Primary  clarifier  effluent  BOD  (mg/1) 

Primary  effluent  suspended  solids  (mg/1) 

Mixed  liquor  suspended  solids  (mg/1) 

Volume  of  air  delivered  (1000  ft^/day) 

Average  mech.  aerator  operating  time  (hrs/day) 
This  should  be  computed  by  multiplying  the 
time  of  operation  of  each  aerator  by  its 
horsepower  and  summing  to  get  the  total. 

This  horsepower-hour  total  should  then  be 
divided  by  the  total  aerator  horsepower. 

Recycle  sludge  volume  (MGD) 

Recycle  sludge  suspended  solids  (mg/1) 

Wasted  sludge  volume  (GPD) 

Average  aeration  basin  dissolved  oxygen  (mg/1) 
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APPENDIX  A 


ACTS  - Subroutine  ACTS  first  calls  PCLAR  which  is  the  primary  clarifier 
subroutine.  Then,  ACTS  reads  in  the  aeration  basin  data  and,  using  the 
output  from  PCLAR  as  input  to  the  basin,  calculates  the  basin  output.  If 
required,  subroutine  ACTNIT  is  called  to  calculate  the  basin  nitrification. 
The  output  from  the  basin  is  then  passed  to  subroutine  ACTSEC  (the  secondary 
clarifier  subroutine)  and  the  systems  final  effluent  is  determined  and 
returned  to  MAIN. 
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ACTNIT  - Subroutine  ACTNIT  takes  the  data  passed  from  ACTS,  determines  the 
aeration  basin  nitrification  and  returns  it  to  ACTS. 

ACTSEC  - Subroutine  ACTSEC  takes  the  data  passed  from  ACTS.  Then,  it  reads 
in  the  secondary  clarifier  data,  calculates  the  plant's  final  effluent  and 
returns  it  to  ACTS. 

PCLAR  - Subroutine  PCLAR  may  be  called  by  either  TFILT  or  ACTS.  PCLAR  takes 
the  plants  influent  data  passed  to  it  by  one  of  these  subroutines.  It  then 
reads  in  the  primary  clarifier  data,  calculates  the  clarifier  output  and 
returns  it  to  the  subroutine  that  called  it. 

■ 
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y G LEVEL  21 


PAIN 


CATE 


75296 


1 3/54/56 


steady  STATE  TREATMENT  PLANT  MODEL 
MODEL  TO  SIMULATE  TRICKLING  FILTER 
OR  ACTIVATED  SLUDGE  TYPE  PLANT 

DIMENSION  TYP(4,2> 

REAL  NH3,N03 
COMMON  FLOW 

OATA  TYP<  1,1)  ,TYP(2,l)  ,TYP(3,1)  ,TYP(4,1  I 

1 /' TR  IC ' » 'KL IN  ' * 'G  FI't'LTER'/t 

2 TYP(l,2)»TYP(2»2)»TYP(3»2)*TYP(4,2) 

3 / ' AC  T I ' , ' VA  TE ' » ' D SL','UDGE'/ 

READ  I 1» 100  II TYP, FLOW, BOD, SS,NH3,N03, TEMP 
WRITE(3,110MTYPI  1,1  TYP), 1 = 1,41 
WRITE! 3,1 15 )FLOW,BOD,SS,NH3,N03, TEMP 
GO  TO  ( 10,201 , ITYP 

10  CALL  TFILT(BOD,SS,NH3,N03, TEMP, £301 
20  CALL  ACTS(B0D,SS,NH3,N03, TEMP, 6301 
30  WRITE! 3, 12C )BOD, SS,NH3,N03 
STOP 

ICO  FORMAT! I l , 6F 8 . 0 ) 

110  FORMAT! 1H1,//'  ESTIMATE  OF  PLANT  PERFORMANCE  AFTER  ', 

1 'INSTALLATION  EXPANS  ION' / 16X , 4A4 , • TYPE  PLANT'// 

2 5X, 'PREDICTED  FUTURE  INFLUENT  CGNC  I T I CNS • / ) 

115  FORMAT! 10X, 'FLOW  R A TE* • , T 3C ,F 8. 2 , 2 X , ' ( MGO ) • / 

1 10X, • BOC* ' »T30»F8.2»2X, • (MG/L ) •/ 

2 10X, • SUSPENDED  SOL  I DS  = ' , T 30 , F 8 . 2 , 2 X , • ! MG/L ) • / 

3 10X, ' AMMON IA='*T30*F8.2»2X, ' (MG/LI •/ 

4 10X,'NITRATE=',T30,F8.2,2X,  • f MG/L ) •/ 

5 10X, 'TEMPERATURE*' ,T30,F8. 2, 2X,’ (DEGREES  C)') 

120  FORMAT! ///5X, 'ESTIMATED  FUTURE  EFFLUENT  CONOITICNS'/ 

2 10X,'B0D=',T30,F8.2,2X,' (MG /LI'/ 

3 10X, 'SUSPENDED  SOL IDS* • , T30 , F8 . 2 , 2X , •( MG/L I • / 

4 IOX, 'AMMONIA*' , T30,F8.2,2X, • (MG/L  I ' / 

5 IOX, 'NITRATE*' ,T30,F 8.2, 2X, • (MG/L I • I 
ENO 


non 


SUBROUTINE  TF I LT ( BOD , SS . NH3 ,N03 , TEMP , ♦ ) 

COMMON  FLOW 

REAL  NH3tNC3,K2C»NN 

DATA  YTFSS  /0.2/ 

CALL  PCLAR ( BOO « SS I 

READ« I, 110JV, A,Rl ,R2,K20tNN 

WRITE(3tI20IVfA»RItR2»K20tNN 

W=B0D*FL0W*8.34 

R=R 1 /FLOW 

F=(1.0«-RI/(  1.0^.1*R»**2 
E=l.0/<  l.0*0.0085*(W/(V*FM**.5> 

E*E* I. 038** (TEMP- 20. > 

IF ( E .GT.  • 95 ) E = « 95 
E*l.-E 

SS=SS*YTFSS*(B0D*8. 34*FL0W-E*B0D*8. 34* ( FLOW+Rl ) I / [ 8. 3** ( FLOWER  I ) ) 

80D*B0D*E 

Q=  l FLOW+-R 1 ) / A 

D*V/A 

I F ( NH3  .GT.  O.JCALL  TF NI T ( NH3 , N03 , TEMP , K20 , NN , D f Q ) 

CALL  FTSEC ( SS »R 1 »R2 ) 

RETURN1 

110  FORMA T ( 6F 8. 0 ) 

120  FORMAT! //1CX, 'TRICKLING  FILTER  CHARACTERISTICS'/ 

1 1 5X » ' TOT AL  FILTER  VOLUME* • , T40 , F 8. 2 . 2X , • ( ACRE-F T ) • , / 

2 15X, 'TOTAL  FILTER  ARE A* • , T40 ,F8. 2 ,2X , • ( ACRES > • » / 

3 15X, 'RECYCLE  TO  FILTER*' »T40,F8.2t2X»' I MGD l • * / 

4 15X, 'RECYCLE  FROM  F I L TER* • , T40 ,F8 .2 ,2 X, • ( MGD I ' / 

5 15X, 'MEDIA  COEFFICIENT  I K ) * • , T40 ,F8.2/ 

6 15X, 'MEDIA  COEFFICIENT  ( N ) * • , T40 , F8.2 ) 

END 


o o o 


LEVEL  21 


MAIN 


CA  IE 


75296 


1 3/54/56 


SUBROUTINE  TFN  I T ( NH 3 , NO 3 , T E , K 2C , NN , n , C ) 
COMMON  FLCM 

REAL  NH3,NC3,K20,K20,NN 
K2U=K20*l.C7**{ TEMP- 20. ) 
E=EXP(-K2C*C/C**NN » 

NO  3=NO  3 + NH3-NH  3*F 

NH3=NH3*E 

RETURN 

ENO 


V G LEVEL  21 


FAIN 


CATC 


752  -*(■ 


13/5A/5t 


SUBROUTINE  TFSEC ( SS , R l , R 2 ) 

COMMON  FLOh 
RE AO ( I* 10 ) A 
WRITE(3,100)A 
FLO*FLOH*Rl-R2 

RF=t l.0-0.82*EXP(-FLC  / ( 2 7 8C . 0* A ) ) t / { l . C- 0. CC2 05*E XP ( -FL C / 
1 (27ft0.0*At)l 

SSXSS*RF 
RETURN 

10  FORMAT(FR.C) 

100  FORHATI//10X, 'SECONDARY  CLARIFIER  CHARACTERISTICS'/ 

1 1 5X  » ' ARE  A= • t T40»F8.2»2X» • (ACRESI * I 

ENO 


n o o 


V G LEVEL  21 


MAIN 


H A TF  = 752™  13/E<./5f 


SUBROUTINE  AC T $ ( ROD , SS , NH 3 , NO  3 , TEMP , * J 

COMMON  FLOW 

REAL  ML  SS ♦ NH  3 * NO  3 

CALL  PCLAR(BOO.SS) 

REACI L,IOO)VOL,MLSS,R,YN 
WRI TEI 3tl2C )VOL,MLSS,R,YN 
PL I=BOO*FLOW/(MLSS*VOL) 

IFIPLI  .GT.  0.*>WRITE(3, 130JPLI 
THETA  = l.O/(YN*PLI  ) 

FLOWR=FLOW+R 

BCD* BOC-ML S S*VOL /! YN* THETA *FLOWR  ) 

SS=MLSS 

IFINH3  .GT.  O.ICALL  AC T N I T ( NH 3 , NOT , F L C W R , VC L , T HE T A , T E M P ) 

CALL  ACTSFCI SS.FLOWR I 
RETURN1 

ICO  FORMA  T ( 4F  8 . C ) 

120  FORMATI //1CX, • ACTIVATEO  SLUDGE  CHARACTERISTICS*/ 

1 15X,* TOTAL  AERATION  BASIN  VOLUME* * » T5C , F 8 . 2 , 2 X , • ( MG  ) * / 

2 15X, 'MIXED  LICUOR  SUSPENDED  SCL I CS = • » T GO » F 8 . 2 » 2 X , • ( M G /l  )•/ 

3 15X,'T0TAL  SLUDGE  REC YCL E = • , T 5C , F B . 2 , 2 X , • ( MGD ) • / 

4 15X  » ' NET  YIELD=* , T50 , F 8 . 2 , 2X  , • (LBS  SCLIDS  PROCLCFO/*/ 

5 T62  » * LB  BOD  UTILIZED)* ) 

130  FORMAT! /5X,  10(**'),'PLI  = ',F5.2/ 

1 15X» • TH IS  PLI  MAY  CAUSE  POORLY  SETTLING  SLUDGE'/ 

2 1 5X , • A ND  INABILITY  TO  MAINTAIN  ESTIMATED  MLSS') 

ENC 


o o o 


V G LEVEL  21 


MAIN 


TATE 


75  295 


1 3/5A/5 f 


10 


22 

2A 

25 


27 


29 


30 
32 

31 
35 


37 


50 

AO 


SUBROUTINE  AC  TN I T ( NH  3 , N03  , F LC  HR  , VC).  , THE  T A , T E MP  ) 

COMMON  FLOW 

REAL  NH3  »NC3»NKS*NKR 

DATA  Y,NKS/0. 05,1.0/ 

R1=FL0WR-FL0W 
TEMPl=VOL/THETA 
NKR=0.33*1. 123** ( TEMP-20. ) 

DT= VOL /FLOWR 

X2 l = Y*F  L0W*NH3  *THET  A/ VOL 
XL  0W  = 0 • 

XUP=2.*X21 
IX  = 0 

S2l=.5*NH3 

SL0W=0. 

SUP=NH3 
IS  = 0 

S22=S21*(NKR*S2l*X21*OT) /( Y*<  NKS+S21 ) ) 

X22=X2l-(NKR*S2l*X21*DT)/(NKS*S2l) 

S17=(S22*FL0WR-NH3*FL0W)/Rl 

X17=X22*FL0WR/R1 

CS2=S? 1-S 1 7 

IF (ABS(DS2) -0.01 >29,29, 22 
I F ( IS— 20>2A,29,29 
IF (DS2  >27,27,25 
SL0W=S2 1 

S21=S2l*0.5*( SUP-SLCW) 

IS=IS+l 
GO  TO  15 
SUP  = S2  1 

S2l=S21-0.5*( SUP- SLOW > 

IS=IS+l 

GO  TO  15 

X l=FL0WR*X2  1 

X2=FLOWR*X22 

X2=(Xl-X2)/X21 

CX2*  TEMP  1-X2 

IF(ABS(CX2>-0.1>50,50,30 

IF  (X2l-0.0C1)A0,A0,32 

t F ( IX- 20) 31, AO, AO 

IF(DX2)37,37,35 

XUP=X2l 

X2l=X2l-0.5*( XUP-XLOW) 

I x = I X ♦ I 
GO  TO  10 
XL0W  = X2  l 

X21=X2l*0.5*( XUP-XLOW) 

I X = I X ♦ I 

GO  TO  10 

NO  3=  NO  3*NH  3- S2  l 

NH3  = S?  I 

CONTINUE 

RE  TURN 

FND 


| 


O O o 


V G LEVEL 


21 


MAIN 


CATE  = 75296 


13/54/56 


SUBROUTINE  ACTSEC ( SS.FLONR) 

COMMON  FLOW 
READ ( i * 100  I A 
WR  I TE  ( 3, 1 10  )A 
A* A*4047. 

SSF*SS*FLOWR*3. 785/24. /A 

SS*4.5*8.6*SSF 

RETURN 

100  FORMAT ( F8 • 0 ) 

110  FORMAT(//IOX, 'SECONDARY  CLARIFIER  CHARACTERISTICS'/ 
l 15X, 'AREA*' ,T40,F 8.2, 2X,» (ACRES  I • ) 

END 


O O O 


V G LEVEL  21 


MAIN 


CATE  = 75296 


SUBROUTINE  PCL AR ( BOD , SS > 

COMMON  FLO* 

READ! 1 • 10 ) A 
WRITE! 3* IOC )A 

RF*(l. 0-0. 82*EXP( -FLOW/ ( 27  80. 0*A > > » / | l. 0-0. C0205*EXP( -FLOW/ 
l (2780.0*A>>> 

SS=SS*RF 

8OD*O.3*BOC+O.7*BO0*RF 

RETURN 

10  FORMAT ( F8 . C ) 

100  FORMAT (//1CX, 'PRIMARY  CLARIFIER  CHARACTERISTICS'/ 

1 15X  t • AREA*'  ,T40,F8.2,2X,' (ACRES)'  ) 

END 


13/5<»/56 
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I 

APPENDIX  B 

Output  From  Treatment  System  Model  Program 

I 

I 

i 

1 

I 

4 

I 
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estimate  of  plant  performance  after  installation  expansion 

TRICKLING  FILTER  TYPE  PLANT 
PREOICTED  FUTURE  INFLUENT  CONDITIONS 


FLOW  RATE* 

14.60 

(MGD) 

BOD* 

148.00 

(MG/L  ) 

SUSPENDED  SOLIDS* 

163.00 

(MG/L) 

AMMONIA* 

0.0 

(MG/L ) 

NITRATE* 

0.0 

(MG/L  ) 

TEMPERATURE* 

15.00 

(DEGREES  C) 

PRIMARY  CLARIFIER 

CHARACTER 

ISTICS 

AREA* 

0.84  (ACRES) 

TRICKLING  FILTER  CHARACTERISTICS 

TOTAL  FILTER  VOLUME* 

4.56 

(ACRE-FT) 

TOTAL  FILTER  AREA* 

1.52 

( ACRES  ) 

RECYCLE  TO  FILTER* 

0.0 

( yGC  ) 

RECYCLE  FROM  FILTER* 

0.0 

( M g 0 ) 

MEDIA  COEFFICIENT  (K)= 

0.23 

MEDIA  COEFFICIENT  (N)= 

2.36 

SECONDARY  CLARIFIER  CHARACTERISTICS 

ARE A=  0 . A 3 (ACRES) 


ESTIMATED  FUTURE  EFFLUENT  CONDITIONS 

BOD*  39.10  (MG/L) 
SUSPENDED  SOLIDS*  6.69  (MG/L) 
AMMONIA*  0.0  (MG/L) 
NITRATE*  0.0  (MG/L) 


ESTI*ATb  OF  PLANT  PERFORMANCE  AFTER  INSTALLATION  EXPANSION 
TRICKLING  FILTER  TYPE  PLANT 

PREDICTED  FUTURE  INFLUENT  CONDITIONS 


FLOW  RATE= 

2.G‘ 

(MOD) 

BOD  = 

250.00 

(MG/L) 

SUSPENDED  SOLI DS= 

500.00 

(MG/L  ) 

AMMON! A = 

20.00 

(MG/L  ) 

Nt  TRATE  = 

0.0 

(MG/L) 

T EMPtR ATU°  E = 

15.00 

(DEGRFES  C) 

PRIMARY  CLARIFIER  CH AR AC TF R I S T I C S 

AR  E A = 

0.50 

( ACRES) 

TRICKLING 

FILTER  CHARACTERISTICS 

TOT  Al 

FILTER  VOL  UM  E = 

5.00 

( ACRE-FT) 

TOTAL 

FILTER  ARF  A = 

C.76 

( ACRES  ) 

RECYCLE  TO  F I L TE R = 

2.o; 

I MGO ) 

TOTAL 

RECYCLE  FROM  FILTE 

2.00 

( MGD ) 

media 

COEFFICIENT  C K > = 

C • 23 

MEDIA 

COEFFICIENT  (N)= 

2.36 

SECONDARY  CLARIFIER  CHARACTERISTICS 

AR  E A=  0.5>  (ACRES) 


ESTIMATED  FUTURE  EFFLUENT  CONDITIONS 

L<On=  7R.7R  (MG/L> 
SUSPENDED  SOLIDS=  15.  >6  (MG/L) 
AMMON  I A=  19. 5R  (MG/L) 
N I TRATE=  2.62  ( MG/L  ) 


APPENDIX  C 

Listing  of  Indicator  Model  Program 

I 

I 

I 

I 

I 

I 

I 

1 

""I 

! 


n o o o o o < > o o o o o o o nan  n o o o o o r>  n 


V L rv/«'L  ?i 


“A  IN 


DATE  = 75296 


98/46/25 


DUX  571 « 

POGFR  P.  S“ITM 


0 1 ME  'AS  l ON  PL  >E SI  4,2) .UNIT (4,4) 

C9mM09  A( 9,6) , PL  ANTI  10) 

OAT  A PLDES(1*2),PLDES(2»2)»PLDES(3»2)»PLDES(4»2) 

1 /MCTI • , 'VATE • , '0  SL ' * * UDGE • /, 

2 “LDESI l, 1 > ,PLDES( 2, 1 » ,PLOES( 3,1 ) ,PLQES(4,l > 

3 /'TRIC'.'KLINS'G  El • * ' L T E R • / 

DA  TA  ’J  J[  T(  1,1  ) ,UNl  T(2,l  ) ,UNI  T ( 3 , 1 ) , UN  I T ( 4 , 1) 

1 /'PRIM*, 'ARY  'SETT*, 'LING'/, 

? UN  I T C 1*2),  UN  (T(2,2I,  'JNI  T ( 3,21,  UN  l T ( 4 , 2 ) 

3 /•  AER' , 'ATIOS 'N  TAS'NKS  '/, 

4 JN  I T ( 1 , 3 ) , UN  I T ( ? , 3 I , UN  I T ( 3 , 3 1 , UN  l T ( 4 , 3 » 

5 /'FINAS'L  CL' ,' ARIF  S' IER  •/ 

DATA  UNI T(  1 ,4) ,UNI T( 2,4 | ,UNI T < 3,4» ,UNIT (4 ,4) 

1 /•  TRlCS 'KL  INS 'G  F l L TER'/ 


S AD I L ,901  ) NOTYP, PLANT 
WRIT'S  3,902) 

WRITE!  3,908)PL ANT 

«■>  1 TE  ( 3,905)  ( PLDFS  ( I ,NPTYP  ) , I = 1 ,4) 
WRITE! 3,903) 

DETERMINE  PROCESS  TYPF 
fin  TO  ( 90 , 80 ) , NP  r YP 

ACTIVATED  SLUDGE  TYPE  PLANT 


r ) READ!  1,906) (All, 11,1  = 1,7) 


PR IMARY  CLARIF  IER 

A ( 1 , 1 ) = TQT AL  TANK  VOLUME  (1000  GAL) 

A ( 2 , l ) = TOTAL  SURTACE  AREA  (FT2) 

A(  3, l ) = TOT  AL  WEIR  LENGTH  (FT) 

A ( 4 , l I =DES I GN  FLOW  RATE  ( M G 0 ) 

A ( 5 , H =DES I GN  DETENTION  PERIOD  (HRS) 

A( 6, L )=DESIGN  SURFACE  SETTLING  RATE  ( GPD/FT2 ) 
A(7,l)=DFSIGN  WEIR  OVERFLOW  RATE  (GPD/ET) 


V.»  I TE  ( 3,904)  ( UNIT!  I , 1 ) , I = 1 ,4) 


IT ( A( 1 , 1 ) .LT. 
I F ( A ( 2 , 1 ) .LT. 
IF  I A ( 3, 1 ) .LT. 

I r ( a c 4 , i » .lt. 
I F I A ( 5 , l ) . L T . 
IF  ( A ( 6 , 1 I .LT. 
ir(A(7,l)  .LT. 


0 . 0 ) WR I TF (3,910) 
0 • 0 I WR  I TE  ( 3 , 91  1 ) 
0.  0 ) WR  I TE  ( 3 , 91  2 ) 
0.0) WRITE (3, 913) 
0.0)WRITE(3,914) 
0. 0 ) WR I TE ( 3,915) 
0 . 0 ) WR I TF ( 3,916) 


READ! 1 , 906 ) ( A ( I ,2)  ,1=1,7) 

AERATION  BASIN 

A(l,2»=T0TAL  TANK  VOLUME  (1000  GAL) 
A( ?,2)=DFS!GN  FLOW  PATE  (mod) 

A( 3,2  I =DE SIGN  DETENTION  PERIOD  (HRS) 


00000010 
00000020 
00000030 
00000040 
00000050 
00000060 
0000007u 
00000050 
000000  JO 

ooooonu 
000001 10 
00QO0120 
00000130 
00000140 
00000150 
00000 160 
00000170 
00000180 
000001  JO 
00000290 
009002  l 0 
00000220 
000002  30 
00000240 
00000250 
00000260 
00 0002 70 
000002*0 
000002  JO 
00000300 
00000310 
00000320 
00000330 
00000340 
00000360 
00000360 
C0000370 
00000380 
00000390 
C0000400 
00000410 
00000420 
C0C004  30 
00000440 
00000450 
C00C0460 
00000470 
000004  80 
OOOC0490 
00000500 
00000510 
CC000520 
000005  30 
00000540 
00000550 
00000560 
00000570 
0C0005H0 


•>  o 


V (1  LEVEL  ?l 


MAIN 


DATE 


75296 


08/66/25 


I 


I 


C 

c 

n 

' * 

C 

r 


C 

C 

C 

c 

C 

r. 

c 

c 

r 

C 


C 

c 

r. 

c 

c 


A(6,2)=DESIGN  ORGANIC  LOADING  (LB  BOO/DAY/IOOO  FT3) 
A(5,?)=DESI GN  *LSS  (MG/L) 

A( 6,2 ) =NUM8ER  OF  AERATORS  OR  COMPRESSORS 
A( 7,2)=TOTAL  AERATOR  OR  COMPRESSOR  HORSEPOWER 


WRITE ( 3*906) (UNI T ( I ,2 ) , I =1 ,6) 

IF ( A ( 1,2)  .LT.  O.OIWRI TE(  1,920) 
I F ( A ( 2 , 2 ) .LT.  0. 0 ) WR  I TE  ( 3 , 92  1 ) 
I F ( A ( 3,21  .LT.  O.OIWRI TF ( 3,922 ) 
ir(A(A,2)  .LT.  O.OIWRI TE( 3,9231 
ir(A(5,2I  .LT.  0. 0 ) WR I TE ( 3 , 926 ) 
I F ( A ( 6 , 2 I .LT.  0. 0 1 WP.  I TE  ( 3 , 925  I 
I F ( A ( 7 , 2 I .LT.  0. 0 1 WR I TE ( 3 , 92  6 I 


R r A 0 ( 1,906  I ( A( I , 3 1 , 1 = 1, RI 


FINAL  CLARIFIER 

A ( 1,31=  TOT  AL  TANK  VOLUME  (1000  GAL  I 
A ( 2 , 3 I =TQT  AL  SURFACE  AREA  (FT2I 
A( 3,3)=TOTAL  WEIR  LFNGTH  (FT) 

A(6,3)=DESIGN  FLOW  RATE  (MGDI 
A ( G , 3 I =OES I GN  DETENTION  PERIOD  (HRS) 

A ( 6, 3 I =DE  SI GN  SURFACE  SETTLING  RATE  (GPD/ET2) 
A ( 7 , 3 ) =DES I GN  WEIR  OVERFLOW  RATE  (GPD/FT) 

A ( M , 3 ) =QES I GN  SOLIDS  LOADING  (LB  SS/DAY/FT2) 


W’R  I TE  ( 3,906  ) (UNIT!  I , 3 > , 1 = 1 ,6  I 
IE ( A ( 1,3)  .LT.  0.0)WRITF(3,910) 
IF ( A I 2 , 3 ) .LT.  0. 0) WRI TF ( 3,91 1 ) 
1 F ( A ( 3 , 3 ) .LT.  0.0  )WR  I TF ( 3,91 2 I 
1 F ( A ( 6 , 3 ) .LT.  0.0)WRITF(3,913) 
I F ( A ( 6 , 3 ) .LT.  O.OIWRI TE(  3,916) 
ir(A(6,3)  .LT.  0.0)WPITE(3,915) 
I F ( A ( 7 , 3 ) .LT.  0.  0 ) WR I TF ( 3 , 91 6 ) 
I F ( A ( 0 , 3 ) .LT.  O.OIWRI TE(  3,91 7) 


GO  TO  103 


TRICKLING  FILTER  TYPE  PLANT 
90  RFAO( 1,906) ( A( I ,1 ) , 1=1 ,7) 


C 

r 

r 

c 

c 

c 

c 


PRIMARY  CLARIFIER 

A(l,ll=TnTAL  TANK  VOLUME  (1000  GAL) 
M2,U*T0TAL  SURFACE  AREA  (FT2) 

A ( 3 , 1 ) =TOT AL  WEIR  LFNGTH  (FT) 

A ( 6 , 1 ) =DE S l GN  FLOW  RATF  ( MGD) 

A ( 5, 1 ) = DE S I GN  DETENTION  “ERIOD  (HRS) 

A(6,l )=DESIGN  SURFACE  SETTLING  RATE  (GPD/ET2) 
A I 7,1) =DE  S I GN  WEIR  OVERFLOW  RATE  (GPO/FT) 


WRITE  13,906 ) ( UNI  T ( ( , l ) , ! = I ,6  ) 


I)  ( A ( 1 , l 
I E ( A ( 2 , 1 
IE ( A( 3,  l 
IF ( A(6,  1 
If  ( A( 6, 1 


.LT.  0.0IWRITF(3,910) 
.LT.  0.0) WR I TE( 3,91 1 ) 
• L T.  0.01  WRITE! 3,912) 
.LT.  O.OIWRI TF( 3,91 3) 
.LT.  O.OIWRI TE( 3,916) 


00000590 
00000600 
00000610 
00000620 
C0C006  30 

0 C000660 
00000650 
00000660 
CC000670 
00000680 
00000690 
C0C00700 
00000710 
00000720 
00000730 
00000760 
00000750 
00000760 
C0CC077O 
00000760 
00000790 
CC000800 
00000810 
00000820 
00000830 
00000860 
CC000850 
00000860 
000008  7 0 
CC0008H0 
C0000890 
00000900 
0CG0091 0 
00000920 
000009  30 
CC000960 
C OU  0095  0 
00000960 
00000970 
00000980 
CC000990 
00001000 
C0001010 
QC0C1020 

oonoio  30 
00001060 
0000 1050 
00001060 
00001070 
00001060 
00001090 
00001 100 
C0O01 1 10 
00001 120 
00001 1 30 
00001 160 
COOOl  1 o 

00001  l 6 u 


O O O O o O O O O Cl  O O o o o o o o o o 


v r.  LEVEL  21 


vain 


DATE  = 75296 


08/46/25 


IF(A(6,1»  .LT.  0.0)WRITE!3,915) 

IF(A(7,1I  .LT.  0.0IWRITEC 3,916) 

C 

flF  AD  ( 1 ,906)  ( A(  I ,2)  , 1=1  ,5) 

TRICKLING  FILTER 
A( 1 * 2 ) =DE  S l GN  FLOW  RATE  (MGD) 

AI 2, 2 )=OESIGN  ORGANIC  LOAD  (LB  BOD/ACRF-FT  DAY) 
A(3,2)= TOTAL  FILTER  VOLUME  (ACRE-FT) 
A(4,2)=T0TAL  FILTER  AREA  I ACRES ) 

A ( 5,  2 ) = TO T AL  RECYCLE  TO  FILTER  (MGD) 

WRITE!  3,904)(UNIT( I ,4)  ,1=1  ,4) 

I F ( A ( 1,2)  .LT.  0.0)WR!TE(3,940) 

I F ( A ( 2 , 2 ) .LT.  O.OIWRITEI 3,941) 

IF (AC  3,2)  .LT.  0.0)WRITF(  3,942  ) 

I F I A ( 4, 2 1 .LT.  0.9)WRITF(3,943) 

I F ( A ( 5 , 2 ) .LT.  0.9)WRI TE( 3,944) 

READ!  1,906) !A(  I ,3)  ,1  = 1,8) 


FINAL  CLARIFIER 

A!  1 , 3 ) = TOT AL  TANK  VOLUME  (1000  GAL) 
A(2,3)=TOTAL  SURFACE  AREA  (FT2) 

A(3,3)=TOTAL  WEI'’  LENGTH  (FT) 

A! 4, 3)=0ESIGN  FLOW  RATE  (MGD) 

A( 5, 3)=0ESIGN  DETENTION  PERIOD  (HRS) 

A ( 6 , 3 ) =DES I GN  SURFACE  SETTLING  RATE  (GPD/FT2) 
A( 7, 3)=. DESIGN  WEIR  OVERFLOW  RATE  ( GPD/FT  ) 

A ( 8, 3 ) =DE S I GN  SOLIDS  LOADING  (LB  SS/DAY/FT2) 


C 

c 


WRITE!  3,904 K UNI T( I , 3) , 1=1  ,4) 

I F ( A ( 1,3) 

.LT. 

O.OIWRI TF! 3,910) 

IF ! A(2,3) 

.LT. 

0 . 0 ) WR I TF ( 3,911) 

IF ( A! 3, 3) 

.LT. 

0. 0 ) WR  1 TE  ( 3 , 91  2 ) 

IF! A(4,3) 

.LT. 

O.OJWRl TF ( 3 , 91 3 ) 

IF! AI5.3) 

.LT. 

0.0)WRITE ( 3,9141 

I F ( A ( 6 , 3 ) 

.LT. 

0.0)WRITE(3,915) 

I F ( A ( 7 , 3 1 

• LT. 

0. 0 ) WR I TF ( 3,916) 

IF! A(8,3) 

.LT. 

O.DtWRI TF (3,917) 

GO  TQ  105 

CALL  ACT!  f. 200) 

105  CALL  TF I L T ( C200 ) 

200  CONTINUE 
STOP 

901  FORMAT! I l , 10A4) 

902  FORMAT! 1 H 1 / / / T4 1 , 50 ( • «" )/ 

1 T4 1 » • * PLANT  DESIGN  AND  OPERATION  CHECK 

2 /T41,50( •*• ) ) 

903  FORMAT! //T20, ‘DESIGN  INFORMATION  NOT  AVAILABLE'/ 

1 T20,32( • * ' ) ) 

904  FORMAT! /T1S,'**'  » 4 A 4 , ' * * ' ) 

905  FORMAT ( T52.4A4, • TYPE  PLANT' I 

906  FORMAT ( 8F 8 . 0 ) 

908  FORMAT! T46 , 10A4 ) 

910  FORMAT ( T20, • TOTAL  VOLUMF') 


00001170 
00001180 
00001190 
00001200 
00001210 
C0001220 
0C001230 
00001240 
C0001250 
00001260 
00001270 
00001280 
00001290 
00001300 
00001310 
00001320 
00001330 
00001340 
COOOl  350 
00001360 
00001370 
COOOl  360 
00001390 
000014 00 
00001410 
00001420 
00001430 
00001440 
00001450 
00001460 
C0001470 
00001480 
00001490 
C0001500 
00001510 
00001520 
00001530 
00001540 
00001550 
C0001560 
00001570 
00001560 
00001590 
CCOC 1600 
00001610 
00001620 
COOO 16  10 
00001640 
00001650 
00001660 
00001670 
00001680 
00001690 
00001700 
00001710 
00001720 
00001730 
00001 740 


V C.  L C VeL  21 


MAIM 


cate 


75296 


08/46/25 


511 

FORMAT! T20  * 

• TOTAL  SURFACE  AREA*  1 

00001750 

712 

TORMAT ( 72  0 » 

• TOTAL  WEIR  LENGTH*  1 

00001760 

913 

rOPM/vx  ( 720, 

• FLOW  RATE  • 1 

00001770 

914 

FORMAT! T20, 

•DETENTION  TIME* 1 

00001780 

915 

FORMAT ( T20, 

•SURFACE  SETTLING  RATE') 

00001793 

916 

FORMAT (T20, 

•WEIR  OVERFLOW  RATE* 1 

CC001830 

917 

FORMAT! T20, 

• SOLIDS  LOADING' I 

0000  1 8 l 0 

92  3 

FORMAT ! T20 , 

•TOTAL  TANK  VOLUME') 

00001820 

921 

FORMAT ( T20, 

' T LOW  RATE*  ) 

00001830 

922 

FORMA  T ( T20, 

•DETENTION  PERIOD' ) 

00001840 

92  3 

TORmat! 720, 

•ORGANIC  LOADING') 

00001850 

92  A 

FO  <m  A7  ( 72  3, 

•MLSS*  I 

00001860 

925 

FORMAT! 720, 

• NUMBER  OF  AERATORS  OR  COMPRESSORS'! 

00001870 

926 

FOR*  A 7 ( 720, 

•TOTAL  AERATOR  OR  COMPRESSOR  HORSEPOWER' 1 

00001880 

74  7 

FORMAT ( T20 , 

•FLOW  RATr* ) 

CC001 890 

94  1 

FORMAT! T20, 

•ORGANIC  LOADING' 1 

00001900 

94? 

FORMAT ( T20, 

• TOTAL  F I L TER  VCLUMF*  I 

00001910 

762 

FORMAT ( T 2 0 , 

• TOTAL  FILTER  AREA'  I 

00001920 

94  4 

FORMAT! T20, 

•FILTr-R  RECYCLE'  ) 

000019)0 

end 

00001940 

n o o o ooooooooor>o*'>  r»  oooor>oooooo  o o o 


V G LFVEt  21 


ACT 


OATE 


75296 


08/46/25 


SUBROUTINE  ACT(*I 
COMMON  A(9,6I .PLANT! 10) 

READ  IN  OPERATIONS  OATA 


RE  AO ( 1 , 901 ) (A!  I ,5)  ,1  = 1 ,8) 


GENERAL  OPERATIONS 
A ( l » 5 ) =PE  AK  FLOW  (INFLUENT)  (MGD) 

A(2,5)=T0TAL  DAILY  FLOW  (INFLUENT)  (MG) 

A ( 3 , 5 ) ■=  BOD  (INFLUENT)  (MG/L) 

A(4,5)=SUSPENDE0  SOLIDS  (INFLUENT)  (MG/LI 
A ( 5 » 5 ) = ROD  (EFFLUENT)  ( MC/l  ) 

A(6,5)=SUSPENDED  SOLIDS  (EFFLUENT)  (MG/L) 

A( 7,5) =AVERAGE  TEMPERATURE  (EFFLUENT)  (DEG  C) 
A(8,5)=D.O.  SATURATION  IN  AERATION  BASIN  (MG/L) 


rtRlTEI 3,909) 
WRITE! 3,910) 

IF  ( A ( 1 , 5 ) .LT. 
I F ( A ( ? , 5 ) .LT. 
IF ( A(  1, 5)  .LT. 
I F ( A ( 4 , 5 ) .LT. 
I F ( A ( r>  , 5 ) .LT. 
I F ( A ( 6 , 5 ) .LT. 
IF ( A( 7,5)  .LT. 
IF ( A(P,5)  .LT. 


0.0)WRITE(3,91l) 
0.0) WRITE (3, 912) 
0. 0) WRI TE ( 3,91 3) 
0 . 0 ) WR I TE (3,914) 
0.0)WRITF(3,915) 
0.0)WRITE(3,916) 
0.0)WRI TF ( 3,91 7) 
0.0) WRITE (3, 918) 


RE  All  ( l ,931  ) ( A ( I ,6)  , 1 = 1 ,9) 


PROCESS  OPERATIONS 
M l ,5) =PRI MARY  EFFLUENT  BOD  (MG/L) 

A(?,6)= PRIMARY  EFFLUENT  SUSPFNOED  SOLIOS  (MG/L) 

A(  3,6)=MIXEO  LIQUOR  SUSPENDED  SOLIOS  (MG/L) 

A ( 4 , 6 ) ^VOLUME  OF  AIR  DELIVERED  (1000  FT3/DAY) 

A ( 5 , 6 ) =HOUR  S Or  OPERATION  OF  MECHANICAL  AERATORS  (HRS/DAYI 
A(6,6)=RETL!RN  SLUDGE  VOLUME  (MGD) 

A ( 7 , 6 ) =R  E TURN  SLUOGE  SUSPENDED  SOLIDS  (MG/L) 

A ( ° , 6 ) =W ASTE  SLUDGE  VOLUME  (MGD) 

A ( 9,6 ) =AVFRAGE  AERATION  RASIN  DISOLVED  OXYGEN  (MG/L) 


WRITE (3,920) 

IF ( A ( 1,6) 

• LT. 

I F ( A ( 2 , 6 ) 

• L T. 

IT  ( A ( 3 , 6 ) 

.LT. 

IT (4(4,6) 

• LT. 

IF  ( A(  5, 6) 

.LT. 

IF ( A ( 6 ,6 ) 

.LT. 

IF ( A( 7,6) 

.LT. 

I T 1 A ( 8,6) 

• L T. 

I F ( A ( 9 , 6 ) 

.LT. 

0. 0) WRI TE ( 3 ,92 l ) 
0.0) WR I TE( 3,922) 
0.0)WRITF(3,923» 
0. 0) WRI TE 13,924) 
0 . 0 ) WR I TE ( 3,925) 
O.OWRJTF  ( 3,926) 
0.0)WRITF( 3,927) 
0. 0 ) WR I TE ( 3,928) 
0.0) WRITE (3,929) 


WR I TE ( 3, 902  ) 

WRITE (3*903)PLANT 
WRITE! 3,904) 


OPERATIONS  AND  OFSIGN  CHECK 


00001950 
00001960 
00001970 
00001980 
00001990 
00002000 
00002010 
00002020 
00002030 
00002040 
00002050 
00002060 
00002070 
00002080 
00002090 
00002 100 
C0002110 
00002120 
00002130 
00002140 
00002150 
00002160 
0C002 1 70 
C00021H0 
00002190 
C00022u0 
C00022  1 0 
00002220 
000022  30 
C0002243 
00002250 
00002260 
00002270 
00002280 
00002290 
C00023 00 
00002  310 
00002  320 
CC002  3 30 
00002  340 
00002350 
00002360 
00002  370 
00002380 
00002  390 
00002400 
00002410 
C0002420 
000024  30 
00002440 
00002450 
00002460 
000024  70 
0C00248G 
000024 90 
00002500 
C000251 0 
000025/C 


J 


o o o o o o no 


DATE  *=  75296 


08/46/25 


| 


V LEVEL  21  ACT 

i 

C PRIMARY  CLARIFIER 

C 

C DETENTION  TIME 

IFIAIl,!)  .LT.  0.  ) GO  TO  101 
XA=All,l)/IA!2,5)*1.0E03/24.) 

XP  = A ( l , 1) /( A( l ,5)*l. OE 03/2  4.  ) 

WRITE! 3 * 9 3 0 ) A ( 5,1  I ,XA, XP 

SURFACE  SETTLING  RATE 

101  I F ( A ( 2 , l ) .LT.  0. I GO  TO  102 
XA=A!2,5»*1.0E06/A!2,l) 

XP  = A I 1.5)*l.0E06/A(2,l  ) 

WR ITF I), 931  )A!6, 1) , XA, XP 

WTIR  OVFRFLOW  RATF 

102  I F ( A ( T , 1 ) .LT.  0. ) GO  TO  103 
XA=A(?,5)*l.OE06/A! 3,1  ) 

XP  = A( 1 ,5)*l.0E06/A( 3,1  ) 

WR  I TC { 1, 932) A ( 7, 1 ) , XA,  XP 

SOLIDS  LOADING 

103  I F ( A ( A , 5 ) .LT.  0.  .OR.  A(2,l»  .LT.  0.)G0  TO  104 
XA=A|4,5)*A(2,5»*R.34/A!2,l) 

NRtHf4»tMAIX*tI”'.34/A(2,  1) 

EFFICIENCY 

104  I F ( A ( 4 , 5 ) .LT.  0.  .OR.  A(2,6)  .LT.  0 . ) GO  TO  105 
EFF  = ( ( At  4 , 5 ) - A ( 2,6)  ) /A  I 4 , 5 ) ) * 1 00. 

WRITE (3*934)EFF 

C AFRATION  PASIN 

C 

C DETENTION  TIME 

105  I F ( A ( 1.2)  .LT.  0.  » GO  TO  106 
I F C A ( 6 , 6 ) .LT.  0.  ) GO  TO  106 
XA=A( 1,2) /( !A!2,5)+A!6,6) ) * l . 0E03 ) *24 . 

XP  = A< 1,2) /( (A(1,5)*A(6,6) ) *1 . 0E03 ) *24  . 

WRITE! 3,940)A( 3»2),XA,XP 

C 

C VOLUME  TR I C LOADING 

106  I F ( A l 1,6)  .LT.  0.  .OR.  All, 2*  .LT.  0.)GO  TO  107 
XA=A!l,6)*A(2,5)*R.34/A!l,2)*7.48l 
XP=A( l ,6) * A I l ,5)*8.34/A( 1,2 )*7.4B1 
wo [ TE  ! 1,94  U XA , XP 

C 

C PROCESS  LOADING  INTENSITY 

IF  I A!  3,6)  .LT.  0. )00  TO  107 

X A = A I 1 ,6) *A ! 2 * 5 > *R  . 34/( A! 3 , 6 ) * A I 1 ,2)*1.0F-03*8.34) 
XP=A(l,6)*A(l,5)*R.34/(A(3,6)*A(l,2)*1.0E-03*8.34) 

WRITE!  3 , 942 ) XA  , XP 

C 

C AVERAGE  SLUDGF  AGE 

107  I F ( A I 3,6)  .LT.  0.  .OR.  All, 2)  .LT.  0.)G0  TO  10R 
I F I A I 8 , 6 ) .LT.  0.  .OR.  AIT, 6)  .LT.  0.)G0  TO  10B 
I F I A I 6 , 5 ) .LT.  0 . ) GO  TO  10R 

XA  = ( A! 3,6) *A! 1 ,2 )*1 .0E-03*fl. 34 1/ ( ( A 18,6 )*1 . E -06* A I 7,6 ) *8 . 34  ) ♦ 
l (IA(?,5)-A!8,6)*l.r-06)*A(6,5)*R.34)) 


00002530 
C0002540 
00002550 
00002560 
CC0025  7 0 
00002580 
00002590 
00002600 
00002610 
00002620 
000026 30 
00002640 
00002650 
00002660 
00002670 
00002680 
C0002690 
00002700 
00002710 
CC002720 
000027  30 
00002740 
00002750 
00002780 
00002780 
00002  790 
00002800 
00002810 
00002R20 
000028  30 
00002840 
00002850 
00002860 
C0002870 
000028H0 
00002890 
00002900 
00002910 
00002920 
00002930 
00002940 
00002950 
00002960 
00002970 
00002980 
C000299C 
00003000 
00003010 
00003C20 
0000  30  30 
00003040 
00003050 
00003060 
00003070 
OCOO  3080 
0000  30  70 
00003100 


I 


i 

i 

I 

1 

I 

I 

i 


LEVEL  ?l  ACT  CATE  = 75296 

wRircn,9<mxA 

c 

C DIFFUSED  AERATION 
C VOLUMrTRIC  AIR  SUPPLY 

105  IF(A(4,6)  .LT.  0.0)  GO  TU  110 
I F I A I 1,2)  .LT.  0. ) GO  TO  109 
XA  = A I 4,6) *109  0.*  7. 4 8052/ 1 & C 1 . ? 1 ♦ 2 <►  - * 6 0 . I 
WRITE!  3,944  ) X A 
C 

C ORGANIC  AIR  SUPPLY 

109  I F ( A ( 1 , 6 ) .LT.  0 . ) GO  TO  11? 

X4  = A( 4,6>*IC00. /( A I 1 ,6)*AI2 ,5 ) *8.34) 

WRITE! 3,945)XA 

GO  TO  112 
C 

C MECHANICAL  AERATION 
C V0LUMFTR1C  POWER  SUPPLY 

in  I F I A ( 7 , 2 ) .LT.  0.0  .OR.  All, 2)  .LT.  0.0)00  TO  112 
XA=A<  7, 2) /A! 1, 2) *7.45052 
WRI TF(3,946)XA 

C 

C DROA  4IC  OXYGEN  SUPPLY 

111  IF ( A ( 9 , 6 ) .LT.  0.  .OR.  All, 6)  .LT.  0.)G0  TO  11? 

I F ( A ( 5 , 6 ) .LT.  0.  .OR.  AIR, 5)  .LT.  O.IGO  TO  112 
ON=?.S*{  | A ( 8, 5)- A I 9, 6)  ) /9. I 7 ) * l . 024** I A ( 7 , 5 ) -20 . ) *0 . 8 
X A = ON*  A l 7 , 2 ) ♦ A I 5, 6 ) / I A I l , 6 ) * A I 2 » 5 ) *8 . 34 ) 

V.R  I TF  I 3 , 94  7 ) X A 

0 

C FFFICIFNCY 

11?  T T I A I 1,6)  .LT.  0.  .OR.  A ( 5 , 5 ) .LT.  O.IGO  Tn  113 
EFr=(All,6)-A(5,5))/A(l,6)*100. 

WP I TF I 3, 9AH ) FFF 

0 SECONDARY  clarifier 

r 

C DETENTION  TIME 

11)  IT  I A ( 1,3)  .LT.  0.  IGO  TO  114 
IF  I Al 6,6)  .LT.  0.  IGO  TO  117 
XA-A(l,3)/(IA(2,5I+A(6»6))*1.0E03)*24. 

X p = A I 1 , 31/1  I A ( l ,5I*A(6,6) ) * 1 .0F03)*24. 

WR  I TE  I 3, 950  A I 5, 3)  , XA,  X° 

C 

C SURf  ACF  SETTLING  RATE 

114  I F | A I 2 , 3 ) .LT.  0.100  TO  115 

XA  = ( (A(2,5)*A(6,6) )*l.0E06) / A I 2 , 3 ) 

X P = I I Al l , 5 ) ♦ A ( 6 , 6 ) l*l.0E06)/AI2»3) 

WRITE!  3 , 95  1 1 A I 6 , 3 » ,XA,XP 

0 

C WEIR  OVERFLOW  RATF 

115  IF ( A I 3, 3)  .IT.  0.  I OH  TO  116 
XA=(A(?*5)+A(6f6) )*1.0F06/A( 3,31 
XP=(A( 1*5)*AI6*6) )*l.0F06/AI3,3) 

W»I  TE I 3,95?) Al 7, 3)  ,XA,XP 

C 

C SOLIDS  LOADING 

116  IF  ( A I 3)  .LT.  O.IGO  TO  117 
IF  I A I 3,61  .LT.  O.IGO  TO  ll« 
XA-AI3,6)*(AI2,5I*A(6,6))*P.34/A(2,3) 


08/46/25 

00903110 
000D112G 
000031 30 
00003140 
COOO  3160 
0000  3 160 
000031  10 
COO  0 3 1 a 0 
0000  3 1 (0 
000O  3200 
00003210 
0000  322  0 
0000  32  30 
0000 J24G 
00003260 
00003260 
000032  7 0 
C00032B0 
00003290 
C0003300 
CCOO  3310 
00003320 
C 0 3 0 3 3 30 
CCOO  3 34  0 
00003350 
0CC0336O 
C00C3370 
000 0 3 3U  0 
00003  390 
OCOO  3400 
000034 1 0 
COOO 3 420 
OCOO  34  30 
000034^0 
OCOO  3460 
00003460 
CCOO  34  70 
0000  34  RO 
COOO  3490 
CCOO  3600 
00003510 
00003620 
CCOO  35  30 
00003640 
0000  3560 
0000  3560 
00003570 
uOO  3 3 6h  0 
COOO  3590 
0000  3600 
0000361 0 
00003620 
0000  36  30 
00003640 
CCOO  36  50 
00003660 
0000  36  79 
0000  36(30 


V G LEVEL  21 


ACT 


DATE 


75296 


03/46/25 


r 

r 


r 

r 


C 


XP  = A ( 3 , 6 ) * I A I l,5)+A(6,6)  )*8.34/A(2,3) 
WRITEI3,953)AI3,3),XA,XP 

FPF 1C IENCY 

117  IF!A(6,51  .LT.  0.  ) GO  TO  113 

LFT=(A( 3, 6)-A(6,5))/A(3, 61*100. 

WRITE!  3*  954  ) EFF 

PI  ANT  300  EFFICIFNCY 

113  1 F ( A I 3 » 5 1 .LT.  0.  .0*.  A(5,5)  .LT.  0.1G0  TC  119 
I=(M3,5)-A(5,5))/A(3,5)*100. 

*n I TE  ( 3 , 960  ) E 

PLANT  SUSPENDED  SOLIDS  EFFICIENCY 

119  IF ( A 1 4 • 5 1 .LT.  0.  .09.  A(6,5)  .LT.  O.IGO  TO  120 

t=(A(4,5)-A(6,5))/A(4,5)*100. 

WR I TE ( 3,961  IE 

120  CONTINUE 
Rc  TURN  l 

901  FORMAT! 10F3. 01 


00003690 
00003700 
CC003710 
00003720 
00003730 
0C003740 
00003750 
00003760 
0CC03770 
00003730 
00003790 
00003800 
00003810 
00003820 
00003830 
000  1 3840 
00003380 
00003860 
0000  3870 
00003R30 
COOO  , >0 


AND  OPERATION  CHECK 


* • 


F i I ° M A T ( T20 , • A VTRAGE 
f DRM  A T ( T20, 'AVERAGE 
FORMAT ( T20 , • AVERAGE 
FORMAT ( T20 , * AVERAGE 


902  FURMATI  1H1///.T41 ,50! •*• ) / 

1 T 4 1 , •*  PLANT  DESIGN 

2 /T41 ,50( •*• > I 
503  FORMAT! T46, 10A4I 

904  FORMAT! //T 30, ‘PARAMETER' ,T 74,' OF  SIGN' , T88 AVER AGE ' ,T106,'PEAK' 
1 T30. •♦*****♦**• ,T74, •♦*****• ,T88, •***♦***', T106, '****• 

909  FORMAT! /T2C,' OPERATIONS  INFORMATION  NOT  AVAILABLE', 

1 /T20,  36!  ' *'  I ) 

910  FORMAT! /T13,'**  GENERAL  OPERATIONS  **'» 

FORMA T( T20, • PFAK  FLOW' I 
FORMAT ! T20, • TOTAL  DAILY  FLOW' ) 

I ORMAT (T20, ' AVERAGE  INFLUENT  BOD') 

INFLUENT  SUSPENDED  SOLIDS') 

EFFLUFNT  BOD* ) 

EFFLUFNT  SUSPENDED  SOLIDS’) 

EFFLUENT  TEMPERATURE') 

FURMAT! T20, 'D.O.  SATURATION  IN  AERATION  BASIN') 

FORMAT I/T18,'**  PROOF SS  OPERATIONS  **') 

FORMAT(T20, 'PRIMARY  EFFLUENT  BOD'  I 
FORMAT! T20, 'PRIMARY  EFFLUENT  SUSPFNOFD  SOLIDS') 

FORMAT ! T20,  • MI XEO  LIQUOR  SUSPENDED  SOLIDS') 

FORMAT ! T20, • VOLUMF  OF  AIR  DEL  I VFREO'  ) 

FORMAT  I T20, 'HOURS  OF  OPERATION  OF  MECHANICAL  AERATORS') 

FORMAT ( T2C, 'RETURN  SLUDGE  VOLUME') 

FORMATIT20, 'RETURN  SLUDGE  SUSPFNDED  SOLIDS') 

FORMAT ( T20, 'WASTF  SLUDGE  VOLUME') 

FORMAT ( T20, • AVTRAGE  AERATION  BASIN  DISOLVED  OXYGEN') 
FORMAT1/T10, 'PRIMARY  CLARIFIER  DETENTION  TIME  (HRS) 

1 ,F 10. 3, T100,F 10. 3) 

FORMAT ( /T 10, • PR IMARY  CLARIFIER 
1 T70,F 10. 3,T85,F10. 3, T100,F10 

FORMAT! /T10, 'PRIMARY  CLARIFIER 
1 T70,F10.3,TB5,F10.3,TIOO,F10.3I 

TORMAT(/T10,'PRIMARY  CLARIFICR  SOLIDS 


911 

912 

913 

914 

915 

916 

917 
913 

920 

921 

922 

923 

924 

925 

926 
92  7 
92  3 

929 

930 

931 
9 32 


SURFACE 
,3  ) 

WEIR  OVERFLOW 


T70.F 10. 3, 


SETTLING  RATE  (GPU/FT2) 


RATE  (GPD/FT ) • 


9 33 


LOADING  RATE  (LB  S5/FT2' 


1 '-DAY)  • , T85,F10. 3, T 1 00, FI  0.3) 


934  FOomaT(/T10, 'PRIMARY 
1 T120.F5.  I ,'  T«) 


CLARIFIER  SOLIDS  HF-ovAL  FFFICIENCY* 


OOOu  39  JO 
00003910 
00003920 
000039  30 
/ 00003940 

/)  OCOO  39  5 0 
0000  3960 
000039  70 
COOO  39«0 
00003990 
000 J40D0 
C00040 1 0 
00004020 
COOU40  30 
0 CO  0404  0 
00004050 
00004060 
C00U407C 
CCCC4080 
000040  > 0 
CC0G4 1 u0 
0C0C41 10 
00004120 
CC004 1 30 
CC004 1 40 
00004 1 50 
C0004 160 
T 850C004  1 70 
00004180 
* 00004190 

CC0C4200 
00004210 
0C004220 
000C42 30 
00004240 
00004250 
CC004260 


V r,  LEVFL  21 


AC  T 


DATE 


75296 


08/46/25 


940  FORMAT-  ( / T l 0 » ' AERATION  BASIN  DETENTION  TIME  (HRS)1# 

1 T70,F10.3,T85,F10.3,TID0,F10.3) 

941  FORMAT! /T10, 'AERATION  BASIN  VOLUMETRIC  LOADING 

1 MLR  BOD/  1000  F T3-D  AY  ) • , TR5  , F 10 . 3 , T 100  , F 10.  3 > 

942  FORMATI/TIO,' PROCFSS  LOADING  INTENSITY  (LB  BOD/LB  MLSS-DAY ) ' , 

I T85.F10. 3,T100,F10. 3) 

943  FORMATI/TIO, 'AVERAGE  SLUDGE  AGE  ( D AYS > • , TB 5 , F 1 0. 3 ) 

944  FORMAT! /TIC,' VOLUMETRIC  AIR  SUPPLY  (CFM/1000  FT  3 ) ' 

1 *T85,F10.3) 

945  FORMAT !/T10, 'ORGANIC  AIR  SUPPLY  IFT3/LB  BOD)', 

I TflS,F10.3) 

946  FORMAT! /T10,' VOLUMETRIC  POWER  SUPPLIED  (HP/1000  FT 3 ) • , 

1 T85.P10.3) 

947  FORMAT! /T10, 'ORGANIC  OXYGEN  SUPPLY  (LB  02/LB  BOD)', 

1 T85.FI0.3) 

94 8 FORMAT ( /T 10, • AERATION  BASIN  BOO  REMOVAL  FF F I C I F NC Y • , T l 20 , F 5. 1 , 

I • r* » 

950  FORMAT! /TIC,' SECONDARY  CLARIFIER  DETENTION  TIME  (HRS)', 

I T70,F10.3,T85,FI0.3,T100,F10.3) 

961  FORMAT! /TIO,' SECONDARY  CLARIFIER  SURFACE  SETTLING  RATE  <GPO/', 

I 'FT2)',T70,F10.3,T85,F10.3,T100,F10.3) 

95?  F0RmaT(/T10,' SECONDARY  CLARIFIER  WEIR  OVERFLOW  RATE  (GPD/FTI', 

I T70.F1 J.3,T85,F10.3,T100,F10.3) 

953  FORMAT! /TIC,' SECONDARY  CLARIFIER  SOLIDS  LOADING  RATE  (LP  SS/FT2', 
1 '-DAY) • , T70,F 10. 3, T85.F10.3, T 100, FI  0.3 ) 

954  FORMAT(/T10,' SECONDARY  CLARIFIER  SOLtDS  REMOVAL  EFFICIENCY', 
l T120,F5.1,'  IT  • I 

960  FOP“AT( /TIC, 'OVERALL  PLANT  BOD  REMOVAL  EFFICIENCY', 
l T l 2 0 , F 5 . 1 , • * • > 

961  FORMAT! /TIO, 'OVERALL  PLANT  SOLIDS  REMOVAL  EFFICIENCY', 

1 T123.F5.1,'  <•) 

END 


C00042  70 
C C 0 042  80 
00004290 
C0004300 
00004  3 I C 
00004320 
00004330 
00004340 
00004350 
00004360 
C0004  370 
C00043H0 
00004390 
C00C4400 
00004410 
00004420 
0 00  044  30 
00004440 
00004460 
00004460 
000044  70 
OCO  044  80 
CC004490 
00QU4-)  30 
00004510 
CC00452  0 
000046  30 
00004540 
00004560 
00004560 
00004570 
CC004580 


V i;  LEVEL  2 I 


TF  !L  T 


DATE 


75296 


08/46/2  5 


I 

I 

I 


C 

c 

r 

C 

r. 

c 

c 

c 

c 

c 

c 

c 


c 

c 

r 

C 

C 

C 

c 

c 

c 


SUBROUTINE  TFILTU) 

COMMON  A(  9,6»  .PLANT!  10) 

READ  IN  OPERATIONS  DATA 

READ! I ,901 ) ! A! 1,5), 1=1,6) 

GENERAL  OPERATIONS 
A ( 1,1) = PF  AK  FLOW  ( INFLUENT)  (MGD) 

A ( 2 , 5 ) =TOT  AL  DAILY  FLOW  UNFLUENT)  ! MG  ) 

A ( 3,5) =BQD  (INFLUENT)  (MG/L) 
A(4,5)=SUSPENDE0  SOLIDS  (INFLUENT)  (MG/L) 
4(5,5) =000  ( EFFLUFNT  ) (MG/L) 

A ( 6 , 5 ) » SUSP END ED  SOLIDS  (EFFLUENT)  (MG/L) 


WRITE*  3,909) 
WRITE!  3,910 
I F ( A ( 1,5)  .LT. 
I E ( A ( 2 , 5 ) .LT. 
I F ( A ( 3,5)  .LT. 
I F ( A ( 4 , 5 ) .LT. 
IF ( A ( 5 , 5 ) .LT. 
I F ( A ( 6 , 5 ) . L T . 


O.OlWRI TF( 3,911) 
0.0)WRITE(3,912) 
0.  0 ) WR  I TE  ( 3 ,91  3 ) 
0.0)WRITE(3,914) 
0.0)WRITF(3,915) 
0.0) WRI TE(3  ,916) 


READ! 1 ,901 ) ( A( 1 ,6  » , 1=1,8) 


PROCESS  OPERATIONS 


A( 1,6) = PRI MARY  EFFLUFNT  BOD  (MG/L) 

A( 2,6)=PRIMARY  EFFLUENT  SUSPENDED  SOLIDS  (MG/L) 
A ( 3,6) =F  I L TER  EFFLUENT  BOD  (MG/L) 

A ( 4 * 6 ) =F  l L TER  EFFLUENT  SUSPENDED  SOLIDS  (MG/L) 
A( 5,6)=RECYCLE  TO  FILTER  ( MGO » 

A( 6,6) =RECYCLF  FROM  FILTER  ( MGD ) 

A ( 7,6 ) ^RECYCLE  FROM  SFCONDARY  CLARIFIER  (MGD) 


WRITE! 3,920) 

IF ( A ( 1 , 6 ) .LT. 
I F f A ( 2 , 6 ) .LT. 
I F ( A ( 3 , 6 ) .LT. 
IF ( A ( 4 , 6 I .LT. 
IF ( A ( 5 , 6 ) .LT. 
IF ( A ( 6 , 6 ) .LT. 
I F ( A ( 7 , 6 ) .LT. 


0.0IWRITF(3,92l) 
O.OIWRI TF ( 3,922) 
0. 0 ) WR I TE ( 3,923) 
0.0)WRITE 13,924) 
0. 0 ) WR  I TE  ( 3 ,92  5 ) 
0. 0 ) WR I TE I 3,926) 
0.0)WRITE (3,927) 


WRITE! 3,902) 

WRITF(3,903)PLANT 
WR I TF (3,904) 

C 

C OPERATIONS  AND  DFSIGN  CHECK 

r 

C PRIMARY  clarifier 

c 

C DETENTION  TIME 

IF ( A( 1 , l ) .LT.  0. )G0  TO  100 
X A - A I 1,1)/IA(2,5)*1. OF  03/24.) 
XP  = A( 1 , 1 »/( A I 1 ,5)*l.0E03/?4. ) 
WRITE! 3,930)4(5,1 ),XA,XP 


C000459C 
00004600 
C0C0461 0 
00004620 
00004630 
C0004640 
C00G46S0 
00004660 
00004670 
00004680 
00004690 
00004700 
C0004710 
00004720 
00004730 
CC004740 
00004760 
00004760 
CC004770 
00004780 
00004790 
CC004800 
00004f>10 
00004820 
C0004B30 
CC004840 
00004850 
C0004860 
CC0C4870 
00004880 
00004890 
C000490C 
00004910 
C0004920 
0C0049  30 
00004940 
00004960 
C 0004  96  0 
00004970 
00004980 
CC004990 
00005000 
00005010 
00005020 
000060  30 
00005040 
00005060 
00005060 
00005070 
00005080 
C0006090 
00005100 
000051 10 
00005120 
CC0C51 30 
00005140 
00005  150 
00005160 


noon  no  no  no  noon  on 


7 0 LEVFL  21 


DATE  = 75296 


08/46/25 


SURFACE  SETTLING  RATE 

100  IF ( A ( 2 * l ) .LT.  0.100  TO  101 
XA=A( 7»5)*l.0E06/A( 2, 1 ) 

X P = A ( 1 ,5)*1.0E06/A(2,1 » 

WRITE!  l,<miAI6,  l)  * X A » XP 

WFIR  OVERFLOW  RATE 

101  |F(A(3,ll  .LT.  0. ) GO  TO  102 
XA=A(2, 51*1. 0E06/A! 3,11 

XP  = A ! I ,51*1 .0E06/A! 3, 1 1 
WR I TE ( 3, 9 32) A! 7,1 ) ,XA,XP 

SOLIDS  LOADING 

102  I F ( A ( 4 , 5 1 .LT.  0.  .OR.  A(2,l)  .LT.  0.)GC  TO  103 
XA=A!4,5)*A(2,5)*R.34/A(2,li 

XP=A!4,5)*A!1 , 5>*8.34/A(2, 1 ) 

WRITE! 3 , 933 1 X A , XP 

EFFICIENCY 

103  I F ( A ( 4 , 5 ) .LT.  0.  .OR.  A(2,6)  .LT.  0.)G0  TO  104 
EFF=I ( A! 4,5)-A!2,6) ) / A ( 4 , 5 1 1*100. 
WRITE!3,934)FFF 

TRICKLING  FILTER 

HYCRAULIC  LOADING 

104  I F ( A ( 4 , 2 ) .LT.  0.)G0  TO  106 
IF(A( 5,6)  .LT.  0. ) GO  TO  106 
XA= ! A! 2, 5) +A( 5,6)  ) / A (4,2  ) 

X P - ( A ( l ,5)  *-A!  5,6)  ) / A ! 4 , 2 ) 

105  XP- ! A!  1, 2)  «-AI  5,2)  ) /A(4,2) 

I F ( A ( l , 2 ) .LT.  0.  .OR.  A ( 5 , 2 ) .LT.  0.)XD=-1. 
WRITE!  3, 9401X0, XA,XP 

ORGANIC  LOADING 

106  I F ( A ( 1 1 6 ) .LT.  0.  .OR.  A(3,2)  .LT.  O.IGO  TO  107 
XA=! A!2,5)*A!  1 , 6 ) * 8. 34 ) /A ( 3 ,2 ) 

XP=( A! 1 , 5 ) * A ! l,6)*8.34)/A( 3,2) 
wRITF!3,94l)A!2,2),XA,XP 

RECIRCULATION  RATIO 

107  IF ( A ( 6 , 6 ) .LT.  O.IGO  TO  108 
XA  = A ( 5 , 6 ) /A  I 2 , 5 ) 

X°=A(5,6)/A!1,5) 

1 OR  XD=A!5,2I/A(l,?» 

I T ! A ! 5 , 2 ) .LT.  0.  .OR.  All, 2)  .LT.  0.)XD  = -1. 
WRITE13, 9421X0, XA,XP 

EFFICIENCY 

109  IF (A! 1,6)  .LT.  0.  .OR.  A|5,5)  .LT.  0.)G0  TO  110 
E=!AI1,6)-A!5,5>)/All, 61*100. 

WRI TE! 3,943)E 

secondary  clarififr 

DETENTION  TI*F 

110  I F ( A ( 1 , 3 ) .LT.  O.IGO  TO  111 


CC0051  70 
C00051H0 
00005190 
00005200 
000052 10 
0000522C 
000052  30 
CC00524C 
00005250 
00005260 
CC005270 
00005280 
00005290 
C0005300 
C0005310 
00005320 
CC005330 
00005340 
00005380 
00005360 
00005370 
00005  380 
00005390 
00GC8400 
00005410 
00005420 
000054  30 
00005440 
00005450 
00005460 
C0005470 
00005480 
C0005490 
0C0C5800 
C00055  10 
0C0U552O 
0C0055  30 
00005540 
C0005560 
00005560 
00005570 
00005  5H  0 
00005590 
00005600 
00005610 
00005670 
00005630 
00005640 
00005650 
00005660 
000056  70 
C00056R0 
00005690 
00005  700 
C0005710 
C0005  770 
00005  7 30 
00005740 


o o 


v <;  LEVFL  21 


TrILT  DATE  = 75296  08/46/25 


IF(A(5,6)  .LT.  0.  .OR.  A(6,6>  .LT.  O.IGO  TO  114 
XA=A(l,3)/((A(2,5|*A(5,6)-A(6,6))*l.OE03)*24. 
XP=A(l,3)/((A(l,5)+A(5,6)-A(6,6))*1.0E03>*24. 
WRITEI 3,950)A(5*3) ,XA,XP 


C 

c 


c 


SURFACE  SETTLING  RATE 

111  IF ( A ( 2 ♦ 3 1 .LT.  O.IGO  TO  112 
XA=(A(2,5)*A(5,6)-A(6,6) ) * l . 0E06/ A ( 2 , 3 ) 

X P = ( A ( l»5)*A(5»6)— A(fc,6)  ) * 1 .0E06/A  ( 2 . 3) 

WR I TE ( 3,951 ) 4(6,3  I ,X A,XP 

WEIR  OVFRF  LOW  RATE 

112  I F ( A ( 3 » 3 I .LT.  O.IGO  TO  113 
XA=(A(2,5)*A(5,6)-A16,6))*l.0E06/A(3,3) 

XP=(  A(  l,5)*A(5,6)-A(6»6))*1.0E06/A(3»3) 
wR  l TE ( 3,95214(7,31 ,XA,XP 

SOLIDS  LOADING 

113  IF (4(4,6)  .LT.  O.IGO  TO  115 

1 F { A ( 8 , 6 ) .LT.  0.  .OR.  A(6,6)  .LT.  O.IGO  TO  114 
I F ( A ( 2 , 3 I .LT.  O.IGO  TO  114 

X4=(4(2,5)*4(5,6)-A(6,6) 1*4(4,61*8.34/4(2, 3 I 
XP-(A(1,5)+A(5,6>-A(6,6)>*A(4,6)*8.34/A(2,3) 
WR(TF( 3,95314(8,3) ,XA,XP 


C 


c 

r 


c 

r. 


c 


CFF IC IENCY 

114  I F ( A ( 6 , 5 I .LT.  O.IGO  TO  115 
E=(A(4, 61-4(6,5) 1/4(4,61*100. 

WR I TE ( 3 , 954  I F 

PLANT  ROD  EFFICIENCY 

115  I F ( A ( 3,51  .LT.  0.  .OR.  4(5,5)  .LT.  O.IGO  TO  116 
C= (4(3,5) -A (5, 511/4(3, 51*1 00. 

WRITE ( 3, 960  IE 

PLANT  SUSPENDED  SOLIDS  EFFICIENCY 

116  I F ( A ( 4 , 5 I .LT.  0.  .OR.  A (6, 5 I .LT.  O.IGO  TO  117 
F=(A(4, 51-4(6, 5)1/4(4,51*100. 

WRITE!  3 , 96 1 I E 


117  CONTINUE 
RE  TURN  l 

901  FORMAT (9F8.0) 

902  FORMAT! 1H1///T41 ,50( '♦•  )/ 

1 Ul  .*<  01  AWT  OF  St 


1 

2 


903 

904 


nr-'*  i | in  i //  / ItlilUl  If 

T 4 1 , • * PLANT  DESIGN  ANO  OPERATION  CHECK 

/T41 ,50( •*• ) I 


* • , 


2 

FORMAT ( T46, 10441 

FORMAT! //T30, 'PARAMETER' ,T74, 'DESIGN* , T88, • AVERAGE* ,T 106, • PE/ 
l T30, '*•*****♦*• ,T74, '•*****' , T 88 , • ♦*♦****• ,T106,'**« 

909  FORMAT ( /TZ0*'0PFR4T I ONS  INFORMATION  NOT  AVAILABLE', 

1 /T20,  36( •••  I I 

910  FORMAT! /T18,'**  GENERAL  OPERATIONS 

911  FORMAT!  T20, 'PEAK  FLOW! 

912  f ORMAT( T20, ' TOTAL  DAILY  TLOH'I 


• *•  I 


912 
91  I 

914 

915 
>16 


FORMAT!  T20  , ' PEAK  FLOW  I 
r O0MAT( T20,  • TOTAL  DAILY  TLOW) 

FOPMAT ( T20, ' AVCRAGE  INFLUENT  HOD* I 
r0RMAT(T20,'AVERAGF  INFLUENT  SUSPENDED  SOLIDS') 
FORMAT! T20,'AVFRAGE  EFFLUENT  BOD' I 
fnoMAT(T20, 'AVERAGE  EFFLUENT  SUSPENDED  SOLIDS') 


00005750 
C0005 760 
00005770 
00005780 
00005790 
00005800 
00005810 
00005820 
00005830 
C0005840 
00005850 
C0005860 
0C005870 
000058H0 
C 0 0 0 5 8 9 0 
CC005900 
00005910 
00005920 
CC0059  30 
00005940 
C0005950 
CC005960 
00005970 
00005980 
C0005990 
00006000 
00006010 
0C0G602G 
000060  40 
00006040 
00006050 
GC006060 
00006070 
0C006080 
C0C06090 
C 00  06 l 90 
C0006  1 10 
CC0C6  120 
00006 1 30 
CC006140 
C0006150 
00006160 
000061 70 
CC006180 
00006190 
C 0006200 
CC0062  10 
00006220 
000062  30 
00006240 
00006250 
00006260 
000062  70 
00006280 
00006290 
CC0C6 300 
00006  4 l 0 
00006  320 


V (.  LEVEL 


TF  I L T 


date 


752-9  6 


08/46/25 


21 


920  FORMAT! /T18,'**  PROCESS  OPERATIONS  **•) 

921  FORMAT!  T20,  'PRIMARY  FFFLlJENT  BOOM 

9 22  FORMAT ( T20* ' PRIMARY  EFFIUFNT  SUSPENDED  SOLIDS*) 

923  FORMAT! T?0,*FILTFR  EFFLUENT  BOD') 

924  FORMAT ( T20* 'FILTER  EFFLUENT  SUSPENDED  SOLIDS') 

925  FORMATIT20, 'RECYCLE  TO  FILTER') 

926  FORMAT! T20, 'RECYCLE  FROM  FILTER') 

927  F0RMAT1T20, 'RECYCLE  FROM  SECONDARY  CLARIFIER') 

930  F0RMAT!/T10,' PRIMARY  CLARIFIFR  DETENTION  TIME  (HRS)', 

1 T70, F 1 0«  3 , TB5  ,F  1 0.  3 » T 1 00  , FI  0 • 3 ) 

931  FORMAT! /T10, ' PRIMARY  CLARIFIER  SURFACE  SFTTLING  RATE  (GPU/FT2)', 
l T70,F10.3,T85,F10.3,T100,F10.3) 

932  FORMAT! /Tl 0,' PRI MARY  CLARIFIER  WEIR  OVERFLOW  RATE  (GPO/FT)', 

1 T70,F10.3,T85,F10.3,T100,F10.3) 

933  FORMATI/TIO, 'PRIMARY  CLARIFIER  SOLIDS  LOADING  RATE  (LB  SS/FT2' 

1 ,'-DAY) • ,T85,F10.3,TIOO,F10.3) 

934  FORMAT! /T10, 'PRIMARY  CLARIFIER  SOLIDS  RE“OVAL  EFFICIENCY', 

1 T120.F5. 1 ,'  *•) 

940  FORMAT! /T10,'  TRICKLING  FILTER  HYDRAULIC  LOADING  ( MGAD  I • , 
l T70,F10.3,T85,F10.3,T100,F10.3) 

941  FORMAT(/T10,' TRICKLING  FILTER  ORGANIC  LOADING  (LB  BOD/ACRE-' 

1 , 'FT-DAYI • ,T70,F10.3,T85,F10.3,T100,F10.3  ) 

94?  FORMAT! /T10, 'FILTER  RECIRCULATION  RATIC  ( R FC I RCUL A T I ON / I NF L UE N I ) ■ 
1 T70,F10.3,T85,F10.3,T100,F10.3) 

943  FORMAT(/T10, 'FILTER  BOD  REMOVAL  EFFICIENCY', 
l T120»F5. 1 , ' *') 

950  FORMATI/TIO, 'SECONDARY  CLARIFIFR  DETENTION  TIME  ( HR S I • . 


00006330 
00006340 
00006350 
00006360 
000063  70 
00006380 
00006390 
00006490 
CC0064  10 
00006420 
000064  30 
C CO  064  4 C 
00006450 
CC  00646  0 
C0006470 
00006480 
00006490 
C0006500 
C00065  10 
OC  006520 
0C0065 30 
00006540 
,000066‘<0 
C0^06S60 
00006570 
C0006580 
CC  006  590 


l T70.F10. 3,T85,F10.3»T1 00 , F l 0 . 3 ) 

951  FORMAT ( /T 10, • SECONDARY  CLARIFIFR  SURFACE  SETTLING  R A I F (GPD/FT2)' 
l T 70,F10.3,T85,F10.3,T100,F10.3) 

952  FORMAT! /TIO,' SECONDARY  CLARIFIER  WEIR  OVERFLOW  RATF  (GPO/FT)', 
l T7U,F10.3,T85,F10.3,T100,F10.3) 

953  FORMAT  (/ T l 0 ,' SECONDARY  CLARIFIER  SOLIDS  LOADING  RATE  (IF  SS/FT?' 
1 '-DAY) • , T70,F10. 3, T8S ,F10.3 , TIOO,F 10. 3 ) 

954  FORMAT! /T10,' SECONDARY  CLARIFIFR  SOLIDS  REMOVAL  EFFICIENCY', 
l Tl20,F5.l,'  *•) 

960  F0RMAT(/T1C, 'OVERALL  PLANT  BOD  REMOVAL  EFFICIENCY', 

1 Tl  20, F5.1*'  V) 

961  FOPMAT ( /TIO,' OVER  ALL  PLANT  SOLIDS  REMOVAL  EFFICIENCY', 

1 Tl20,F5.l,'  *•) 

END 


C 0 0066  00 
, 00  0066  I 0 
CCC06620 
000066  30 
00006640 
, C0006650 
00006660 
000066  ro 

CC0066HQ 
COO  066  90 
00006  700 
C0006  7 10 
CC006  720 
00006  7 30 


APPENDIX  D 

Output  From  Indicator  Model  Program 


* plant  ncsi-'A  a.'.n  lpcp.UIC-j  Chl 

FL3T  CRL.  C.\ 

TPlCKLINb  FILTL-l  T YP  r PLANT 


